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FOREWORD 


The  report  "Effluent  Mixing  Zone  Studies"  is  one  of  a  series  of 
technical  documents  prepared  for  the  Grand  River  Basin  Water 
Management  Study.  The  project  described  herein  was  undertaken 
through  the  Grand  River  Study  Team  at  the  request  of  the  Grand  River 
Implementation  Committee. 

The  material  contained  in  these  reports  is  primarily  technical 
support  information  and,  in  itself,  does  not  necessarily  constitute 
policy  or  management  practices.  Interpretation  and  evaluation  of 
the  data  and  findings,  in  most  cases,  cannot  be  based  solely  on  this 
one  report  hut  should  be  analyzed  in  light  of  other  reports  produced 
within  the  comprehensive  framework  of  the  overall  study.  Questions 
with  respect  to  the  contents  of  this  report  should  be  directed  to 
O.G.  Weatherbe,  P.Eng.,  Supervisor,  River  Systems  Assessment  Unit, 
Water  Resources  Branch,  135  St.  Clair  Avenue  West,  Toronto,  Ontario, 
M4V  1P5. 
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EXECUTIVE  SUMMARY 


One  of  the  concerns  of  the  Grand  River  Basin  Water  Managenent  Study 
is  the  need  to  define  the  toxic  effect  areas  of  chlorine  and  ammonia 
from  wastewater  effluent.  This  is  accomplished  by  the  delineation 
of  a  "mixing  zone"  or  "limited  use  zone"  in  a  portion  of  the 
cross-section  of  river  below  the  outfall,  where  concentrations  of 
the  pollutants  may  not  comply  with  the  provincial  water  quality 
objectives  (PWQO).  The  remaining  portion  of  the  cross-section  of 
the  river,  where  the  water  quality  is  suitable  for  maintaining 
aquatic  life,  is  termed  the  "zone  of  passage".  Water  management 
policies  of  the  Ministry  of  the  Environment  (MOE)  allow  the  use  of 
mixing  zones  provided  terns  and  conditions  governing  the  size  and 
water  quality  within  the  mixing  zone  are  met.  These  provisions  are 
stipulated  in  the  MDE  publication  "Water  Management  (MOE,  1978)". 

A  two-dimensional  mathematical  model  was  used  to  determine  limited 
use  zones  and  zones  of  passage  attainable  in  individual  mixing  zones 
of  receiving  streams  under  present  and  future  loading  conditions  for 
the  various  management  options  considered  in  the  long-range  planning 
study.  The  model  incorporates  longitudinal  convection,  lateral 
dispersion  and  first-order  decay  of  a  nonconservative  pollutant.   In 
this  study,  five  wastewater  treatment  plants  were  considered: 
Waterloo,  Kitchener,  Gait,  Rrantford,  and  Guelph.  Effluents  from 
the  first  four  plants  are  discharged  into  the  Grand  River;  the 
Guelph  plant  discharges  its  effluent  into  the  Speed  River.  The 
two-dimensional  model  was  not  applied  to  the  Brantford  case  since 
the  effluent  was  observed  to  spread  across  the  entire  width  of  the 
Grand  River  in  a  short  distance  below  the  outfall.  For  the  other 
four  cases,  the  reductions  in  effluent  concentrations  required  to 
obtain  desired  zones  of  passage  were  calculated  for  specific  design 
conditions  (t  g.  effluent  flow  rates  and  concentrations  of  chlorine 
and  ammonia,  various  river  flow  rates  and  a  lateral  boundary  of  the 
limited  use  zone  equal  to  40%  of  river  flow  rate).  For  the 
Rrantford  case,  concentrations  of  un-ionized  ammonia  and  chlorine 
were  predicted  using  the  instantaneous  completely  mixed  assumption. 
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Results  of  the  mixing  zone  study  indicated  that  within  the  study 
areas  of  the  Grand  River  Basin,  for  the  design  cases  considered,  a 
zone  of  passage  having  a  lateral  boundary  consisting  of  60%  of  the 
total  river  flow  is  attainable  with  regard  to  chlorine  toxicity. 
Only  after  the  2001  medium  projected  population  level  at  Guelph, 
under  low  flow  winter  conditions,  would  the  PWQO  for  chlorine  be  in 
noncompliance  in  the  Speed  River. 

In  order  to  achieve  the  PWQO  for  un-ionized  ammonia,  tertiary 
treatment  is  necessary  at  the  Waterloo  wastewater  treatment  plant. 
With  the  implementation  of  tertiary  treatment  at  Kitchener,  the  PW00 
for  un-ionized  ammonia  can  be  achieved  until  2031  (summer  low 
flow);  however,  with  advanced  treatment  (activated  carbon 
absorption),  no  violations  of  the  PWQO  would  occur.  Current 
tertiary  treatment  at  the  Guelph  plant  is  not  sufficient  for 
maintaining  the  necessary  zone  of  passage  in  the  Speed  River,  based 
on  the  PWQO  for  un-ionized  ammonia.  At  Gait,  no  additional 
treatment  over  the  conventional  activated  sludge  is  necessary  to 
maintain  satisfactory  levels  of  un-ionized  ammonia  in  the  river. 

In  the  case  of  the  Grand  River  below  the  Brantford  wastewater 
treatment  plant,  results  based  on  the  instantaneous  completely  mixed 
assumption  indicate  that  residual  chlorine  would  have  a  minimal 
effect  (extending  approximately  20  m  downstream  of  the  outfall), 
whereas  un-ionized  ammonia  concentrations  may  not  be  in  compliance 
with  the  PWQO  for  downstream  distances  of  124  m  to  2239  m  under 
different  conditions  of  streamflow  and  temperature. 
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I.   INTRODUCTION 

1.1  General 

Wastewater  effluents  discharged  into  receiving  streams  and 
rivers  often  contain  substances  such  as  chlorine  and  ammonia 
which  are  potentially  toxic  to  aquatic  biota.  Other  water  uses 
such  as  recreation,  agriculture  and  withdrawals  for  drinking 
water  supply  can  also  be  affected  by  toxic  levels  of  such 
pollutants.  In  order  to  minimize  the  potentially  toxic  effects 
of  such  pollutants,  effluents  are  generally  discharged  through 
bank  outfalls  so  as  to  create  a  plume  that  hugs  a  shoreline  of 
the  river. 

Due  to  the  geometric  and  hydraulic  characteristics  of  an  outfall 
and  of  a  receiving  watercourse,  transverse  and  longitudinal 
variations  of  concentrations  of  the  toxic  pollutants  occur  as 
the  effluent  moves  downstream  from  the  point  of  discharge.  The 
zone  of  a  river  between  the  outfall  and  the  nearest  cross 
section  where  the  instream  mixed  concentration  attains 
uniformity  is  termed  the  "mixing  zone". 

In  order  to  provide  a  favourable  habitat  for  the  biota  within  a 
mixing  zone,  it  is  desirable  to  maintain  a  specific  portion  of 
the  cross  section  of  the  receiving  stream  wherein  the 
concentrations  of  the  pollutants  comply  with  a  specified  water 
quality  objective  (C  ).  The  portion  of  the  stream  in  which 
the  water  quality  is  considered  suitable  for  the  maintenance  of 
aquatic  life  is  termed  the  "Zone  of  Passage"  (ZOP).  The 
remaining  portion  of  the  river  where  the  pollutant  concentration 
may  not  comply  with  the  specified  objective  is  termed  the 
"Limited  Use  Zone"  (LUZ).  The  salient  features  of  a  mixing  zone 
are  depicted  in  Figure  1  (Gowda,  1980a  and  1980b).  The  policies 
and  guidelines  describing  the  use  of  mixing  zones  is  presented 
in  the  publication  "Water  Management"  (MOE,  1978). 
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Figure  1.      Salient  features  of  a  Mixing   Zone   in  a   River. 

(Note:      Cross-sections    1-4   refer  to   Figure   3) 


One  of  the  concerns  of  the  Grand  River  Basin  Water  Management 
Study  (GRBWMS)  is  the  maintenance  of  zones  of  passages  in  river 
channels  below  the  major  wastewater  treatment  plants  in  the 
megalopolis  region  of  the  Grand  River  Basin.  Of  particular 
concern  are  the  potential  problems  of  chlorine  and  ammonia 
toxicity  to  the  aquatic  biota  which  inhabit  the  rivers  below  the 
treatment  plants. 

This  study  describes  the  application  of  a  two-dimensional  steady 
state  mathematical  model  which  incorporates  advection,  lateral 
dispersion  and  decay  of  a  pollutant  to  predict  instream 
concentration  distributions  of  chlorine  and  ammonia  for  the 
various  management  options  being  considered  in  the  long-range 
river  basin  planning  study.  The  results  presented  in  this 
report  should  be  considered  as  indicative  only  of  general, 
relative  trends  in  water  quality  changes  resulting  from  various 
management  options.  Given  the  data  limitations  outlined  within 
the  report,  detailed  field  and  modelling  studies  at  appropriate 
implementation  stages  are  recommended. 

1.2  Objectives 

The  main  goal  of  this  study  is  to  assess  the  ability  of  the 
various  management  alternatives  to  provide  a  zone  of  passage 
with  respect  to  the  provincial  surface  water  quality  objectives 
for  chlorine  and  ammonia.  The  specific  objectives  of  the  study 
are  as  fol  lows: 

a)  To  determine  if  a  specified  ZOP  for  aquatic  biota  is 
attainable  with  respect  to  toxicity  of  total  residual 
chlorine  and  of  ammonia,  under  present  and  future  loading 
conditions,  by  maintaining  the  present  concentrations  of 
the  pollutants  in  the  effluents. 

b)  To  determine  if  the  ZOP  is  attainable  under  present  and 
future  loading  conditions  for  the  various  management 
options  (producing  various  effluent  pollutant 
concentrations)  considered  in  the  long-range  planning  study. 


c)  To  estimate  the  allowable  effluent  concentration  of  each 
pollutant  so  as  to  achieve  the  required  ZOP. 

d)  To  estimate  the  maximum  longitudinal  boundary  at  which  the 
specified  water  quality  objective  for  each  pollutant  is 
achieved. 

1.3  Pollutants  Considered  and  Their  Toxicity 

Current  wastewater  treatment  practices  for  mechanical  plants 
within  the  Province  of  Ontario  require  effluent  disinfection, 
generally  through  chlorination,  before  an  effluent  can  be 
discharged  into  a  receiving  water  body.  Residual  chlorine 
present  in  the  chlorinated  effluents  is  potentially  toxic  to 
aquatic  biota.  The  toxicity  component  of  residual  chlorine  is 
usually  expressed  in  terms  of  Total  Residual  Chlorine  (TRC), 
which  is  given  by  the  sum  of  free  (H0C1  and  0C1")  and  combined 
(chloramines)  forms.  Ammonia  is  commonly  present  in  municipal 
wastewater  effluents  in  both  ionized  (NH.)  and  un-ionized 
(NHL)  forms,  the  fraction  of  each  being  dependent  on  the  pH 
and  temperature  of  the  receiving  water.  The  un-ionized  portion 
of  total  ammonia  is  potentially  toxic  to  aquatic  biota. 

In  order  to  ensure  a  satisfactory  level  of  water  quality  in 
surface  water  bodies,  the  Ministry  of  the  Environment  (1978)  has 
published  a  set  of  objectives  termed  "Provincial  Water  Quality 
Objectives"  (PWQO).   In  general,  the  instream  concentration  of  a 
pollutant  discharged  into  a  receiving  water  body  should  comply 
with  the  PWQO  for  that  pollutant  to  ensure  protection  of  aquatic 
life,  recreation,  and  other  uses.  When  a  limited  use  zone  is 
utilized  for  dilution  of  the  effluent,  the  specified  instream 
water  quality  objective,  C       must  be  met  at  a  given  lateral 
boundary  of  the  LUZ,  in  order  to  ensure  a  zone  of  passage  for 
fish  and  other  biota.  The  PWQO  for  total  residual  chlorine  and 
un-ionized  ammonia  are   0.002  mg/L  and  0.02  mg/L,  respectively 
(maximum  daily  permissible  concentration). 
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The  Ministry's  Surface  Water  Policy  1  (1978)  states  that  in 
areas  where  water  quality  is  better  than  the  PWQO,  the  water 
quality  must  be  maintained  at  or  above  the  objectives.  Policy  2 
states  that  for  areas  where  water  quality  does  not  currently 
meet  the  objectives,  no  further  degradation  may  be  incurred  and 
all  practical  measures  should  be  taken  to  upgrade  water  quality. 

Therefore,  the  intended  goal  of  any  management  scenario  is  to 
achieve  the  PWQO.  In  view  of  the  fact  that  these  objectives  are 
not  always  achievable,  particularly  for  surface  waters  currently 
in  non-compliance  with  the  PWQO  (as  outlined  in  Policy  2), 
alternative  management  plans  which  could  achieve  an  improvement 
in  surface  water  quality  or  at  the  very  least  do  not  further 
degrade  surface  water  quality  need  to  be  examined. 

1.4  Water  Pollution  Control  Plants 

This  study  considers  individual  mixing  zones  in  rivers  below 
four  major  water  pollution  control  plants  (WPCP)  -  three  on  the 
Grand  River,  viz.,  Waterloo,  Kitchener  and  Gait,  and  the  Guelph 
WPCP  on  the  Speed  River  (Figure  2).  The  present  treatment  level 
and  the  design  hydraulic  capacity  of  each  plant  is  presented  in 
Table  1. 

TABLE  1:  PRESENT  TREATMENT  LEVEL  AND  CAPACITY  FOR  FIVE  WPCP'S  WITHIN 
THE  GRAND  RIVER  BASIN 


WPCP  TREATMENT  LEVEL*        HYDRAULIC  CAPACITY 

MGD        m3/sec 


Waterloo  1  10.0  0.44 

Kitchener  1  27.0  1.18 

Guelph  2  12.5  0.55 

Gait  1  8.5  0.37 

Brantford  1  18.5  0.81 

*  Treatment  Level  1  is  conventional  activated  sludge  only;  and  Levpl 
2,  known  as  tertiary  treatment,  is  a  conventional  activated  sludge 
plant  with  rotating  biological  discs  for  nitrification  and  sand 
filters  for  effluent  polishing.  Effluent  chlorination  is 
practiced  at  each  WPCP. 
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figure  2.       Locations  of  the  Water  Pollution  Control  Plants 
under  study  within  the  Grand    River  Basin. 


The  two  dimensional  nixing  zone  model  utilized  in  this  study  is 
designed  for  a  shore-attached  effluent  plume  below  a  bank 
outfall  in  a  river,  and  has  been  applied  to  the  locations  below 
Waterloo,  Gait  and  Guelph.  The  actual  outfall  for  the  Kitchener 
wastewater  treatment  plant  is  located  in  the  middle  of  the 
river.  Field  studies  carried  out  at  this  location  during  1975 
and  1979  have  shown  that  the  released  effluent  moves  towards  the 
right  bank  (looking  upstream),  encounters  a  riffle,  then  crosses 
the  river  to  the  left  bank  at  a  45  angle,  thus  creating  a 
potentially  toxic  barrier.  The  effluent  then  tends  to  hug  the 
left  bank  as  it  continues  to  move  downstream.  Since  the  mixing 
zone  model  requires  a  bank  outfall,  the  point  of  discharge  has 
been  assumed  to  be  located  at  the  left  shoreline,  below  the 
riffle  (about  300  m  downstream  of  the  present  location)  for 
modelling  purposes. 

In  the  case  of  the  Brantford  Wastewater  Treatment  Plant,  the 
effluent  mixes  completely  across  the  river  in  a  short  distance 
below  the  bank  outfall.  Therefore,  the  two  dimensional  model 
(applicable  to  shore-attached  effluent  plumes)  was  not  applied 
to  the  Brantford  case.  However,  instream  concentrations  of 
un-ionized  ammonia  and  chlorine  have  been  predicted  based  on  the 
instantaneous  complete  mixing  assumption  and  are  presented  in 
this  report. 

The  model  only  considers  the  effects  of  pollutant  concentrations 
in  the  effluent  of  the  wastewater  treatment  plants  under  study. 
Background  concentrations  of  pollutants  were  not  incorporated 
into  the  model  predictions. 
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II  MODELLING  ASPECTS 

2.1  Stream  Tube  Model 

2.1.1  Description  of  Model: 

The  basic  concept  of  the  stream  tube  model,  developed  by 
Yotsukura  and  Cobb  (1972),  is  the  utilization  of  partial 
cumulative  discharge,  q,  at  a  given  cross-section  of  a  river  as 
the  independent  variable  in  the  analysis.  Hence  the  river  is 
divided  into  a  number  of  vertical  strips  termed  "stream  tubes", 
such  that  the  discharge  within  each  strip  is  the  same.  This 
transformation  conveniently  accounts  for  the  solute  mass  balance 
and  adequately  approximates  natural  stream  conditions  resulting 
from  variable  widths  and  velocities  (Yotsukura  ?i   Cobb,  1972). 
Cross-sectional  concentration  distributions,  c(x,  q),  predicted 
by  the  streamtube  model  will  therefore  be  functions  of  q  rather 
than  the  normally  used  transverse  distance,  y.  The  details  of 
the  derivation  of  the  streamtube  model  have  been  presented  by 
Yotsukura  (1977),  Yotsukura  and  Sayre  (1976)  and  Yotsukura  and 
Cobb  (1972).  Modifications  of  the  stream  tube  model  equations 
to  account  for  longitudinal  variabilities  in  decay  rate 
coefficients  and  hydraulic  parameters  of  river  channels  are 
presented  by  Gowda  (1980b).  The  derivations  of  these  analytical 
expressions  are  beyond  the  scope  of  this  report,  but  are 
available  in  the  above  citations.  However,  the  important 
aspects  of  the  model  are  briefly  described  in  the  following 
sections.  All  computational  procedures  are  applicable  to 
shore-attached  plumes  in  shallow  rivers  resulting  from  effluent 
discharge  at  bank  outfalls. 


2.1.2  Limited  Use  Zone  Boundary 

In  order  to  maintain  an  acceptable  zone  of  passage,  the  lateral 
boundary  of  a  LUZ  is  limited  to  some  specified  fraction  of  the 
total  river  flow.  The  partial  cumulative  discharge  between  the 
outfall  bank  and  the  accepted  LUZ  lateral  boundary  is  denoted  by 
q.  .  The  specified  instream  water  quality  objective,  Cs, 
must  be  met  at  the  lateral  boundary,  q.  ,  at  all  downstream 
distances  below  the  outfall.  Instream  concentrations  of  a 
particular  pollutant  for  a  given  q  =  q.  increase  to  a  maximum 
at  some  value  of  x  and  then  decrease  as  shown  in  Figure  3.  The 
point  at  which  the  concentration  reaches  a  maximum  value  is 
termed  the  "critical  point".  The  concentration  at  this  critical 
point  is  termed  the  "critical  concentration",  C,  ,  and  the 
distance  downstream  of  the  outfall  to  the  critical  point  is 
termed  the  "critical  distance",  X.  .  Methods  of  determining 
X.  and  C,  are  described  in  other  publications  (Gowda,  1980a 
and  1980b).  The  co-ordinates  of  a  critical  point  for  a  given 
q,  are  used  to  calculate  allowable  effluent  concentrations  as 
explained  in  the  following  sub-section.  The  longitudinal 
boundary  of  a  LUZ  occurs  at  the  first  downstream  point  along  the 
effluent  shoreline  where  the  water  quality  objective,  C  is 
met. 

2.1.3  Allowable  Effluent  Concentration 


Once  the  critical  concentration,  C,  for  a  given  q.  is  known, 
the  allowable  effluent  concentration,  C  .,  can  be  calculated 
according  to  the  following  equation  (Gowda,  1980a): 


c    =      c    c 

LeA      Le  S 


cL  (1) 


C  =  concentration  of  pollutant  in  effluent 
e 

C  =  specified  instream  concentration  criterion 

C,  =  critical  (maximum  instream)  concentration  for  a  given  q, 
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When  C  .  is  less  than  C  ,  a  reduction  in  the  effluent 
eA  e' 

concentration,  C  ,  is  required.  When  C  .  is  greater  than 

C  ,  the  stream  is  capable  of  assimilating  an  effluent 

concentration  higher  than  C  . 
3        e 

2.2  Computational  Aspects 

A  computer  model  MIXAPPLN,  written  in  FORTRAN  language,  has  been 
developed  using  the  stream  tube  model  equations  to  perform 
design  case  computations.  A  documentation  for  the  program  taken 
from  Gowda  (1980b)  is  provided  in  Appendix  A.  A  brief 
description  of  the  program  is  provided  below. 

Data  on  channel  hydraulics  and  pollutant  decay  rates  obtained 
from  field  surveys,  along  with  design  case  options,  are  input 
into  the  model  to  produce  a  matrix  of  concentrations  for 
cumulative  values  of  q  at  various  distances  below  the  outfall. 
Through  an  iterative  procedure,  the  program  calculates  allowable 
effluent  concentration  and  maximum  longitudinal  spread  (i.e., 
longitudinal  boundary  of  LUZ)  for  a  given  set  of  input  data. 

Concentrations  at  any  distance,  x,  downstream  of  the  outfall  are 
calculated  according  to  the  following  relation: 

C(x,  q)  =  (0.5)(C  )(RF)(SUM)/(tt0)1/2        (2) 


C    =  cross  sectional  average  concentration  in  river  just 

a 

below  outfall. 


RF   =  decay  function,  which  accounts  for  reach-dependency 
of  the  decay  rate  including  temperature  effects. 

6      =   B1  (X)  ,  where 
(B) 
B'  =  nondimensional  coefficient  (accounting  for 

variations  in  channel  geometry  and  transverse 
dispersion  characteristics),  determined  from 
field  data  and  model  calibration  study. 
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X  =  downstream  distance  from  outfall  (m). 
B  =  width  of  channel  at  distance  X  (m). 
SUM  =  a  term  to  account  for  the  effects  of  lateral 

dispersion  and  sidewall  reflection  of  pollutant  on 
C(x,  q) 

Figure  4  depicts  graphically  the  iterative  process  utilized  by 

the  program  to  calculate  C  ..  The  concentration  distributions 

predicted  by  Equation  2  are  utilized  to  search  for  the  transect 

at  which  the  concentration  is  maximum  for  a  given  q,  .  Using 

this  concentration,  C   .,.,  and  the  associated  downstream 

cnt 

distance,  X   .. ,  along  with  the  weighted  mean  values  of  width 
cnt 

and  velocity  as  first  estimates,  new  values  of  C.  and  X.  are 
calculated.  If  the  absolute  percentage  difference  (100  (X.  - 
X   4.)/X,  ),  is  greater  than  an  arbitrarily  chosen 
percentage,  the  procedure  is  repeated  until  the  difference  is 
within  the  specified  limit.  For  the  purpose  of  this  study,  the 
acceptable  difference  was  chosen  to  be  5%.  The  values  obtained 
in  the  final  iteration  step  are  then  printed  out  as  the  critical 
concentration,  C. ,  and  critical  distance,  X.  ,  for  a  given 
q.  .  Knowing  C,  and  Cs,  the  allowable  effluent 
concentration,  C  .,  is  calculated  according  to  Equation  1. 


NOTE:  In  the  program  listing  (Appendix  A),  taken  from  Water 

Resources  Paper  #14  (Gowda  1980b),  the  dimensionless  term  6 
has  been  modified  by  replacing  B  with  BW,  where  BW  denotes  the 
moving  average  value  of  channel  width  for  the  distance,  X, 
downstream  of  the  outfall.  The  coefficient,  R',  was  replaced 
with  a  dimensionless  coefficient,  B",  where  B"  «  B'(BW/B).  A 
comparison  of  the  predictions  of  the  two  computer  program 
versions  for  the  case  of  total  ammonia  in  the  Grand  River 
below  Waterloo  indicated  differences  in  C  .  of  up  to  0.4 
mg/L  (as  N)  and  in  maximum  longitudinal  spreads  of  up  to  300 
meters.  Predicted  values  for  longitudinal  spread  given  in 
Water  Resources  Paper  #14  are  consistently  lower  than  those 
presented  in  this  report.  On  the  other  hand,  predicted 
concentrations  in  the  former  publication  are  higher  than  those 
in  this  report. 
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When  the  effluent  is  discharged  from  a  bank  outfall,  the  maximum 
concentrations  at  various  transects  within  the  mixing  zone  occur 
along  the  outfall  shoreline  (q  =  0).  In  order  to  determine  the 
downstream  distance  along  the  shoreline  where  the  instream  water 
quality  objective  is  met,  the  program  searches  for  the  transect 
where  the  predicted  shoreline  concentration  is  just  above  C-. 
The  program  then  iterates  in  a  similar  manner  as  outlined  above 
to  determine  the  corresponding  distance,  X   ,  where  C<.  is 
met. 

In  addition,  the  model  also  performs  a  similar  iterative 
procedure  to  determine  the  downstream  shoreline  distance 

XsceA'  where  tne  ^stream  water  quality  is  met  when  the 

effluent  concentration  is  C  .. 

eA 

2.3  Scale-up  of  Model  Parameters 

In  order  to  assess  present  and  future  water  quality  within  the 
Grand  River  Basin  with  regard  to  chlorine  and  ammonia  toxicity, 
various  management  options  have  been  examined  (see  Chapter 
III).  These  options  include  treatment  plant  expansion,  flow 
augmentation,  and  seasonal  variabilities  in  streamflow  and 
temperature.  Changes  in  design  temperature  and  design 
streamflow  would  require  a  scale-up  of  model  parameters  (width, 
depth,  velocity,  etc.)  to  account  for  the  variabilities 
associated  with  the  various  options. 

The  decay  rates  of  total  residual  chlorine  and  ammonia  are 
dependent  on  temperature  and  are   adjusted  in  the  MIXAPPLN 
program  according  to  the  modified  Van't  Hoff-Arrhenius 
relationship  given  by 

KT2  =  KT1  0    Ï2"Tl  <4) 
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KT1  ,  KT?  =  decay  rate  coefficients  at  temperatures  T, 
and  T«  °C,  respectively;  and  $    =  temperature  correction 
factor.  Values  of  6    for  TRC  and  ammonia  were  set  equal  to  1.03 
and  1.106,  respectively  (Gowda,  1978;  Eckenfelder,  1970). 


Changes  in  streamflow  rates  result  in  corresponding  changes  in 
channel  width,  depth  and  velocity.  The  computer  model  utilizes 
the  Leopold-Maddock  equations  to  acount  for  the  changes  in 
hydraulic  parameters  (Eckenfelder,  1970).  For  example,  a  width 
R(2)  resulting  from  a  streamflow  Q(2)  is  calculated  by 

B(2)  =  B(l)(  Q(2)PMR  (5) 

U(i)j 

0(1)  =  measured  flow 

R(l)  ■  measured  width  corresponding  to  Q(l) 

BPWR  =  empirically  derived  constant 

A  more  detailed  description  of  scale-up  procedures  is  provided 
in  other  publications  (Gowda,  1980a  and  1980b). 
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Ill  EVALUATION  OF  MANAGEMENT  OPTIONS 

3.1  Options  Model  led 

The  various  management  plans  under  consideration  for  the  Grand 
River  Basin  are  fully  described  in  Smith  (1981)  and  Will  son  et 
al  (1981).  The  various  options  and  the  design  parameters 
associated  with  each  option  are  presented  below  (subsections 
3.1.1  -  3.1.5).  These  values  are  used  as  inputs  to  the  mixing 
zone  model,  along  with  instream  process  parameters  (see 
subsection  3.1.6). 

3.1.1  Effluent  Flow  (QEFL) 

Present  and  future  water  quality  depends  on  changes  in  effluent 
loading  rates  relating  to  projected  population  increases  in  the 
major  urban  centres.  The  effects  of  projected  increases  in 
populations  and  resultant  pollutant  loads  for  20  and  50  year 
planning  periods  beginning  with  1981  (i.e.  years  2001  and  2031) 
are  incorporated  into  the  model  in  terms  of  effluent  flow  rates 
(QEFL).  The  projected  population  for  the  year  2001  is  a  medium 
density,  while  both  medium  and  high  density  populations  were 
projected  for  the  year  2031  (to  account  for  possible  changes  in 
future  growth).  The  effluent  flows  utilized  in  this  study  are 
summarized  in  Table  2. 

TABLE  2:  EFFLUENT  FLOWS  (QEFL) 


WATERLOO 
cfs    m^/s 

KITCHENER 
cfs    m^/s 

GUELPH 

GALT 

BRANT 
cfs 

FORD 

YEAR 

cfs 

mJ/s 

cfs 

mJ/s 

rrvVs 

1979 

10.86 

0.31 

23.70 

0.67 

18.20 

0.52 

9.09 

0.26 

17.35 

0.49 

200 1M* 

21.22 

0.60 

3i   28 

1.11 

27.92 

0.79 

17.43 

0.49 

:i„55 

0.69 

2031 M 

30.62 

0.87 

56.26 

1.59 

50.70 

1.43 

25.02 

0.71 

37.37 

1.06 

2031  H** 

35.21 

1.00 

64.86 

1.84 

66.18 

1.87 

28.77 

0.79 

46.67 

1.32 

*  M:  Medium  density  population 
**  H:  High  density  population 
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3.1.2  Upstream  Flows  (QRUP) 


The  7-consecuti ve-day  low  flow  with  a  20-year  return  period 

(7Q20)  under  natural  unregulated  discharge  conditions  is 

approximately  1.1  m  /sec  (40  cfs)  in  the  Grand  River  at 

Bridgeport.  There  are  two  existing  reservoirs  on  the  Upper 

Grand  River,  namely  Lake  Belwood  and  Lake  Conestoga.  The  river 

3 
discharge  is  currently  augmented  by  approximately  8.5  m  /s 

(300  cfs)  from  the  two  reservoirs  from  spring  until  fall.  The 

3 
river  discharge  is  approximately  1.1  m  /sec  (40  cfs)  during 

winter  since  no  augmentation  occurs  due  to  the  low  reservoir 

levels. 


Six  upstream  flow  rates  in  the  Grand  River  at  Bridgeport  are 
considered  in  the  design  cases,  based  on  the  following 
management  conditions: 

Design  Flow 

1:       1.1   m3/sec   (40  cfs)  winter  flow   (7Q20). 

3 
2:       4.2  m  /sec   (150  cfs)         higher  winter  flow   resulting   from 

warm  wet  winter  or  the  use  of 

reservoirs  during  spring-fall 

strictly  for  maintaining  flood 

control  rather  than  flow 

augmentation. 

3 

3:   8.5  m'  /sec  (300  cfs)    the  present  target  of  350  cfs 

3 
4:   11.3  m  /sec  (400  cfs)   summer  low  flow  at  Doon. 

3 
5:   14.2  m  /sec  (500  cfs)   incorporates  the  construction  of 

the  West  Montrose  reservoir  to 

maintain  water  quality  in  future 

through  flow  augmentation. 
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17.0  m  /sec  (600  cfs)   incorporates  the  building  of  both 

the  West  Montrose  and  the  St. 
Jacobs  Reservoirs  for  flow 
augmentation. 


The  Speed  River  currently  has  one  reservoir  located  above  the 
City  of  Guelph.  Three  river  flows  upstream  of  Guelph  are 
considered  in  the  design  cases: 

Design  Flow 

3 
1:   0.5  m  /sec  (17  cfs)     unregulated  natural  flow  at  Guelph 

(7Q20). 

3 
2:   1.8  m  /sec  (64  cfs)     present  flow  with  augmentation  from 

the  Guelph  reservoir. 

3 
3:   2.5  m  /sec  (88  cfs)     augmentation  from  the  proposed 

reservoir  at  Everton  or  a  50% 

reduction  in  the  permanent  flood 

pool  at  the  Guelph  reservoir. 


When  determining  upstream  flow  rates  for  Kitchener  and  Gait, 
increases  in  the  upstream  WPCP  effluent  flow  rates,  due  to 
projected  increases  in  populations,  must  be  considered.  Also, 
since  the  Gait  WPCP  is  located  below  the  confluence  of  the  Speed 
and  the  Grand  Rivers,  upstream  flow  rates  for  the  Gait  WPCP 
include  a  combination  of  both  river  flows. 

Table  3  presents  the  design  upstream  flow  (QRUP)  utilized  in  the 
modelling  studies  for  each  WPCP. 
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TABLE  3.   UPSTREAM  FLOWS  (QRUP)  UTILIZED  FOR  THE  DESIGN  CASES 


YEAR 

DESIGN 
OPTION 

WATERLOO 

3/      « 
m  /s    cfs 

KITCHENER 

3/         r 

n  /s    cfs 

GUELPH 

GALT 

m3/s 

cfs 

3, 

m  /s 

cfs 

1979 

1 

1.1 

40 

1.4 

50.9 

0.5 

17 

3.2 

114.7 

2 

4.2 

150 

4.6 

160.9 

1.8 

64 

6.4 

224.7 

3 

8.4 

300 

8.8 

310.9 

2.5 

88 

11.9 

421.7 

4 

11.2 

400 

11.6 

410.9 

14.8 

521.7 

5 

14.0 

500 

14.5 

510.9 

18.3 

645.7 

6 

16.8 

600 

17.3 

610.9 

21.1 

745.7 

2001M* 

1 

1.1 

40 

1.7 

61.2 

0.5 

17 

4.4 

154.0 

2 

4.2 

150 

4.8 

171.2 

1.8 

64 

7.5 

264.0 

3 

8.4 

300 

9.1 

321.2 

2.5 

88 

13.0 

461.0 

4 

11.2 

400 

11.9 

421.2 

15.9 

561.0 

5 

14.0 

500 

14.8 

521.2 

19.4 

685.0 

6 

16.8 

600 

17.6 

621.2 

22.2 

785.0 

2031M* 

1 

1.1 

40 

2.0 

70.6 

0.5 

17 

5.9 

206.9 

2 

4.2 

150 

5.1 

180.6 

1.8 

64 

9.0 

316.9 

3 

8.4 

300 

9.4 

330.6 

2.5 

88 

14.5 

513.9 

4 

11.2 

400 

12.2 

430.6 

17.4 

613.9 

5 

14.0 

500 

15.0 

530.6 

20.9 

737.9 

6 

16.8 

600 

17.8 

630.6 

23.7 

837.9 

2031H* 

1 

1.1 

40 

2.1 

75.2 

0.5 

17 

6.7 

237.4 

2 

4.2 

150 

5.2 

185.2 

1.8 

64 

9.8 

347.4 

3 

8.4 

300 

9.5 

335.2 

2.5 

88 

15.4 

544.4 

4 

11.2 

400 

12.3 

435.2 

18.2 

644.4 

5 

14.0 

500 

15.2 

535.2 

21.7 

768.4 

6 

16.8 

600 

18.0 

635.2 

24.6 

868.4 

M  and  H  denote  medium  and  high  density  population  levels,  respectively, 
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3.1.3  Temperature  and  pH 

As  stated  previously  (Section  1.3),  total  ammonia  is  present  in 
wastewater  effuents  in  ionized  (NH.)  and  un-ionized  (NhL) 
forms,  the  latter  form  being  potentially  toxic  to  fish.  The 
fraction  of  each  form  is  dependent  on  pH  and  temperature 
according  to  the  equation  (Willingham  1976): 


C 
u 


1  +  10pKa  "  pH 


pKa  =  0.09018  +  2729.92  6) 


273.2  +  T° 


C    =  concentration  of  total  ammonia 

C    =  concentration  of  un-ionized  ammonia 
u 

pH   =  pH  of  receiving  water 

pKa  =  -log  K  ,  where  K,  =  dissociation  constant  for  total 
a        a 

ammonia 
T°   =  temperature  of  receiving  water,  °C 

pH  values  of  7.8  and  8.2,  based  on  observed  field  data,  were 
utilized  for  the  study.  Temperature  values  used  in  the  model 
range  from  5°C  to  25  C  in  5  degree  steps.  Lower  values  of 
5°  and  10°C  are  considered  representative  of  winter 
conditions,  while  the  range  15  -  25  is  representative  of 
late  spring,  summer  and  early  fall  temperature  variation 
experienced  in  the  basin. 

3.1.4  Treatment  Levels 

All  the  wastewater  treatment  plants  within  the  study  area  of  the 
basin,  with  the  exception  of  the  Guelph  plant,  are  currently 
operating  as  conventional  activated  sludge  plants.  In  order  to 
assess  the  effects  of  various  treatment  levels  on  instream 
concentrations  of  chlorine  and  ammonia  within  the  mixing  zones, 
three  treatment  levels  are  examined  in  the  study: 
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1.  Conventional  Activated  Sludge  (CAS):  Primary  and  secondary 
treatment  to  remove  carbonaceous  BOD;  effluent  disinfection 
to  destroy  micro-organisms  of  public  health  concern. 

2.  Tertiary:  CAS  plant  with  nitrification  to  reduce 
nitrogenous  BOD,  sand  filtration  for  effluent  polishing; 
and  disinfection  of  effluent. 

3.  Advanced:  Tertiary  with  activated  carbon  adsorption  (ACA) 
and  effluent  disinfection.  ACA  reduces  carbonaceous  and 
nitrogenous  BOP,  and  also  tends  to  remove  small  amounts  of 
ammonia. 

Table  4  presents  typical  effluent  concentrations  of  total 
residual  chlorine  and  total  ammonia  attainable  by  the  three 
levels  of  treatment. 

TABLE  4:  EFFLUENT  CHARACTERISTICS  FOR  VARIOUS  WASTEWATER  TREATMENT 
LEVELS 


TREATMENT  LEVEL 


TOTAL  RESIDUAL 
CHLORINE 


TOTAL  AMMONIA 


Conventional  Activated  Sludge 

Tertiary 

Advanced 


500  ug/L 
500  ug/L 
500  ug/L 


12.5  mg/L 
1.5  mg/L 
1.0  mg/L 


3.1.5  Limited  Use  Zone  Boundary 

The  publication  "Water  Management"  (MOE,  1978)  states  that: 

"The  mixing  zone  dimensions  will  be  kept  as  small  as 
possible  while  ensuring  that  the  provincial  water  quality 
objectives  are  met  at  the  boundary."  It  is  also  stated 
that  "a  mixing  zone  shall  not  be  allowed  to  create  a 
barrier  to  the  migration  of  fish  and  aquatic  life." 
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Therefore,  the  specified  inst ream  water  quality  objective,  C  , 
must  be  met  at  an  acceptable  lateral  boundary  of  the  LUZ,  q,  , 
to  ensure  an  adequate  zone  of  passage  for  aquatic  biota. 


In  order  to  examine  the  sensitivity  of  the  lateral  boundary  of 
LUZ  on  the  allowable  effluent  concentration  and  maximum 
longitudinal  spread,  model  runs  were  made  with  q.  ranging  from 
20%  to  40%  of  the  total  river  flow,  Q.  Based  on  preliminary 
analyses  and  discussion  with  Grand  River  Water  Quality 
Sub-Committee  members,  a  LUZ  of  40%  of  the  total  river  flow  was 
considered  appropriate  for  the  Grand  River  Basin.  Consequently, 
the  remaining  60%  of  the  total  river  flow  represents  the  zone  of 
passage.  (Note:  All  model  runs  include  predictions  of 
allowable  effluent  concentrations  and  maximum  longitudinal 
spreads  with  p.  =  0.2,  0.3  and  0.4,  where  p.  =  q./Q). 

3.1.6  Instream  Process  Parameters 

The  instream  process  parameters  required  as  model  inputs  are: 
(i)  streamflow,  channel  width,  mean  depth  and  mean  velocity  at 
known  distances  in  each  river  stretch  below  the  outfalls,  as 
well  as  the  empirical  constants  of  the  Leopold-Maddock 
equations;  (ii)  the  dimensionless  coefficient,  B';  and  (iii) 
first-order  decay  rate  coefficients  of  TRC  and  total  ammonia. 
The  hydraulic  parameters  including  the  coefficient,  B',  were 
obtained  from  field  surveys  carried  out  at  each  of  the  study 
locations  during  1975.  The  decay  rates  of  TRC  were  based  on  the 
values  determined  from  field  data  analyses  in  other  streams  and 
rivers  (Gowda  1978  and  1980a).  The  ammonia  decay  rates  were 
obtained  from  the  results  of  previous  waste  assimilation 
studies.  As  stated  earlier,  the  instream  process  parameters 
were  scaled  up  to  account  for  changes  in  streamflow  and 
temperature  values.  Instream  process  parameters  utilized  for 
this  study  can  be  fi/jnd  in  the  Appendix. 
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The  general  procedures  involved  in  field  data  collection, 
parameter  estimation  and  scale-up  of  various  parameters  were 
described  in  another  publication  (Gowda  1980b),  which  should  be 
consulted  for  details. 

3.2  Modelling  Results 

Figures  5  through  14  illustrate  the  results  of  the  model  runs 
for  total  residual  chlorine  and  un-ionized  ammonia  for  the 
mixing  zones  below  the  Waterloo,  Kitchener,  Guelph  and  Gait 
wastewater  treatment  plants.  As  stated  previously,  the  stream 
tube  model  was  not  applied  to  the  Brantford  case.  Results  of 
runs  with  temperatures  of  5  C  and  25  C  are  presented  on  the 
figures  to  show  the  effects  of  two  temperature  extremes. 
Although  two  pH  options  were  utilized  for  ammonia,  the 
predictions  with  pH  =  8.2  are  represented  on  the  figures,  since 
the  concentration  of  un-ionized  ammonia  is  greater  at  a  higher 
pH,  when  all  other  parameters  are  held  constant. 

A  temperature  of  5  C  is  considered  representative  of  winter 
conditions;  therefore,  plots  of  runs  at  5°C  are  presented  with 
the  corresponding  winter  discharges  on  the  Grand  River  of  1.1 
m3/sec  (40  cfs)  and  4.2  m3/sec  (150  cfs)  at  Bridgeport. 
Plots  of  predictions  with  a  temperature  of  25°C  and  upstream 
flows  at  Bridgeport  of  8.4  m3/sec  (300  cfs)  to  16.8  m3/sec 
(600  cfs)  are  considered  to  be  representative  of  summer 
conditions.  For  the  Speed  River,  the  combination  of  5°C  and 

0.5  m  /sec  (17  cfs)  represent  winter  conditions,  and  25°C 

3  3 

and  1.8  m  /sec  (64  cfs)  and  25  m  /sec  (88  cfs)  represent 

summer  conditions.  Since  present  and  future  water  quality  is 

dependent  on  changes  in  effluent  loading  rates,  the  total 

pollutant  load  given  by  the  product  0  .C  (effluent  flow 

rate  x  effluent  concentration)  expressed  in  m i  1 1  i g ..-.-•  s  per 

second,  is  plotted  against  instream  critical  concentration, 

C.  ,  in  Figures  5  to  14.  The  figures  show  increasing  load 

through  time  at  constant  flow  (x-axis)  and  the  resultant 

critical  concentration  (maximum  at  LUZ  boundary)  achieved  in  the 

stream  (y-axis).  Several  design  upstream  flow  conditions  are 

shown  for  each  case. 
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In  order  to  maintain  a  zone  of  passage  which  is  the  portion  of 
the  channel  containing  60%  of  the  total  river  flow,  the  critical 
concentration,  C.  ,  of  the  pollutant  must  not  exceed  the 
specified  instream  water  quality  guideline,  C-,  at  a  LUZ 
boundary  of  40%  of  the  total  river  flow,  as  previously  stated. 
Therefore,  portions  of  the  curves  extending  above  C-   are 
indicative  of  conditions  where  a  60%  zone  of  passage  is  not 
attainable. 


The  following  subsections  describe  the  results  of  model 

predictions  at  each  location.  For  those  conditions  where  the 

PWQO  objectives  are  not  met,  allowable  effluent  concentrations 

which  would  result  in  compliance  with  the  PWQO  have  been 

calculated.  The  results  are  summarized  in  Tables  5,  6  and  7  for 

the  Waterloo,  Kitchener  and  Guelph  cases,  respectively. 

Downstream  shoreline  distances  (X   .),  where  the  PWQO  for 

un-ionized  ammonia  is  met  have  also  been  calculated;  the  results 

are  presented  in  Table  8.  Generally,  the  smaller  values  of 

X   .  are  associated  with  higher  temperature  and  lower 

streamflow  rates.  The  values  of  X   .  where  the  PWQO  for 

sceA 

total  residual  chlorine  is  met,  were  generally  small  (in 
comparison  with  those  for  un-ionized  ammonia)  and  are  not 
presented  in  this  report. 

In  Figures  5  to  14,  the  instream  critical  concentrations  of  TRC 
for  a  40%  lateral  boundary  of  LUZ  are  seen  to  increase  with  an 
increase  in  streamflow,  whereas  the  critical  concentrations  of 
ammonia  decrease  with  an  increase  in  the  streamflow.  An 
examination  of  the  various  TRC  modelling  runs  indicates  that  the 
critical  concentration  of  TRC  increases  with  streamflow  up  to 
some  value  and  then  declines.  The  reasons  for  such  predictions 
are  not  clear;  however,  different  combinations  of  hydraulic 
parameters  and  decay  rates  have  been  shown  to  result  in  such 
predictions  (Howda  1978;  1980a). 
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3.2.1  Waterloo 

Figure  5  illustrates  the  results  of  the  model  runs  for  total 
residual  chlorine  for  the  mixing  zone  in  the  Grand  River  below 
the  Waterloo  WPCP.  The  curves  do  not  extend  above  the  PWQO  of  2 
ug/L  for  TRC,  indicating  compliance  with  respect  to  chlorine 
toxicity  during  both  winter  and  summer  conditions. 

Figure  6  shows  the  plots  for  un-ionized  ammonia  with  the 
Waterloo  WPCP  operating  at  the  CAS  level  of  treatment.  The  PWQO 
of  0.02  ug/L  for  un-ionized  ammonia  is  exceeded  under  virtually 
all  conditions  of  loading  and  temperature.  The  addition  of 
tertiary  treatment  (in-plant  nitrification)  at  the  Waterloo  WPCP 
would  result  in  un-ionized  ammonia  concentrations  for  both 
winter  and  summer  conditions  being  within  the  PWQO  (Fig.  7). 

For  conditions  where  the  PWQO  for  un-ionized  ammonia  are  not  in 
compliance,  the  allowable  effluent  total  ammonia  concentrations, 
C  .,  necessary  to  limit  the  lateral  boundary  of  the  LUZ  to 
40%,  were  calculated  according  to  Equation  1.  The  results, 
presented  in  Table  5,  show  that  except  for  upstream  flows 
greater  than  14.2  m  /s  (500  cfs)  under  the  1979  population 
conditions,  the  values  of  allowable  effluent  total  ammonia 
concentrations  are  generally  smaller  than  12.5  mg/L  (the  nominal 
total  ammonia  in  a  CAS  treatment  plant).  These  results  indicate 
that  tertiary  treatment  (nitrification)  is  required  at  the 
Waterloo  WPCP. 

The  values  of  X,  .  for  total  ammonia,  presented  in  Table  8, 
sceA 

show  that  the  shoreline  distances  at  which  the  PWQO  is  met  would 
be  significantly  smaller  with  tertiary  treatment  at  the  Waterloo 
WPCP,  in  comparison  to  those  attainable  with  the  CAS  treatment 
level. 
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TABLE  5:   ALLOWABLE  EFFLUENT  CONCENTRATIONS*  OF  TOTAL  AMMONIA  AT  THE  WATERLOO 
WPCP  TO  MAINTAIN  A  40%  LIMITED  USE  ZONE 


Allowable  Effluent  Concentration 
Upstream  Flow     Temperature  to  Meet  PWQO  Targets 

m3/s     cfs        °C        Year  mg/L  (as  N) 


1.1  40         5         1979  4.39 

2001M  2.69 

2031M  2.14 

2031H  1.98 

4.2  150         5         2001M  7.30 

2031M  5.33 

2031H  4.75 

8.5      300        25         1979  7.67 

2001M  4.05 

2031M  2.88 

2031H  2.54 

11.3     400        25         1979  9.97 

2001M  5.22 

2031M  3.69 

2031H  3.24 

14.2      500        25         1979  12.24 

2001M  6.39 

2031M  5.22 

203 1H  3.94 

17.0      600        25         2001M  7.43 

2031M  5.22 

2031H  4.57 


The  effluent  concentration  of  total  ammonia  at  Waterloo  is 
presently  12.5  mg/L.  Allowable  effluent  concentrations  are 
those  required  (i.e.  less  than  12.5  mg/L)  to  achieve  PWQO 
targets. 
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3.2.2  Kitchener 

The  modelling  results  for  total  residual  chlorine  indicate  that 
a  zone  of  passage  which  is  60%  of  the  total  river  flow  in  the 
Grand  River  downstream  of  the  Kitchener  wastewater  treatment 
plant  is  achievable  for  all  design  case  population  levels, 
upstream  flows  and  temperatures  (Figure  8). 

Figure  9  shows  noncompliance  of  the  PWQO  for  un-ionized  ammonia 

for  1979  winter  and  summer  conditions,  with  Kitchener  operating 

at  conventional  treatment.  When  tertiary  treatment  at  the 

Kitchener  WPCP  is  simulated,  the  PWQO  are  in  noncompliance  only 

3 
during  the  summer  conditions  at  8.4  m  /sec  (300  cfs)  after 

approximately  2001  and  by  2031H  for  11.2  m3  sec  (400  cfs),  as 

seen  in  Figure  10.  With  ACA  treatment  at  Kitchener,  the  results 

presented  in  Figure  11  show  that  the  PWQO  for  un-ionized  ammonia 

is  met  under  all  the  design  cases. 

Table  6  summarizes  the  allowable  effluent  concentrations  of 
total  ammonia  necessary  to  achieve  the  PWQO  for  ammonia;  the 
values  range  from  1.28  to  7.11  mg/L.  The  results  indicate  the 
need  for  inplant  nitrification  at  the  Kitchener  WPCP  under  the 
1979  design  condition. 

For  total  ammonia,  the  range  of  X  .  is  11  to  11,650  meters 

and  greater  with  CAS  treatment,  650  to  14,100  meters  with 

tertiary  treatment  and  325  to  6,750  meters  with  ACA  treatment 
(Table  8). 
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TABLE  fi.  ALLOWABLE  EFFLUENT  CONCENTRATIONS  OF  TOTAL  AMMONIA  AT  THE  KITCHENER  WPCP 
TO  MAINTAIN  A  40%  LIMITED  USE  ZONE 


Allowable  Effluent 
Treatment  Present  Effluent  Concentration  to  Meet 

Level     Concentration    Upstream  Flow  Temperature  Year      PWQO  Target 

nw/sec    cfs      °C  mg/L  (as  N) 


Conventional    12.5  mg/L      1.1      40 
Activated      total 
Sludge        ammonia 


4.2     150 


8.5     300      25 


11.3     400      25 


14.2     500      25 


17.0     600      25 


Tertiary       1.5  mg/L      8.5     300      25 

total 
ammoni  a 


1979 

2.95 

2001M 

2.38 

2031M 

2.08 

2031H 

1.98 

1979 

7.11 

2001M 

4.87 

2031M 

3.82 

2031H 

3.48 

1979 

3.96 

2001M 

2.56 

2031M 

1.90 

2031H 

1.69 

1979 

5.01 

2001M 

3.18 

2031M 

2.34 

2031H 

2.07 

1979 

6.03 

2001M 

3.80 

2031M 

2.76 

2031H 

2.44 

1979 

7.04 

2001M 

4.41 

2031M 

3.19 

2031H 

2.81 

2031M 

1.43 

2031H 

1.28 
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3.2.3  Guelph 

There  are  no  problems  in  attaining  a  60%  zone  of  passage  with 
regard  to  total  residual  chlorine  in  the  Speed  River  below  the 
Guelph  wastewater  treatment  plant  (Figure  12)  under  summer  high 
flow  conditions.  Only  after  the  2001M  projected  population 
level  at  Guelph,  under  low  flow  winter  conditions,  are  the  PWQO 
violated. 

With  the  addition  of  tertiary  treatment  (nitrification)  at 
Guelph,  the  PWQO  for  ammonia  is  violated  in  1979  under  summer 
conditions  with  the  current  population  and  by  2001M,  for  winter 
conditions  (Figure  13). 

The  X   .  values  are  in  the  range  of  4025-8200  meters  with 
sceA  3 

regard  to  ammonia  under  tertiary  treatment  conditions  (see  Table 
8). 

The  allowable  effluent  concentration  of  total  ammonia  and 
chlorine  under  various  design  conditions,  have  been  presented  in 
Table  7.  In  the  case  of  total  ammonia,  the  values  of  C  .  are 
lower  than  1.5  mg/L  indicating  a  need  for  further  reduction  in 
effluent  ammonia  concentrations  (through  advanced  treatment 
methods).  In  order  to  meet  the  PWQO  for  residual  chlorine, 
dechlorination  of  effluent  by  the  year  2031  is  also  indicated  by 
the  results  (for  chlorine)  presented  in  Table  7. 
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TABLE  7.  ALLOWABLE  EFFLUENT  CONCENTRATIONS  OF  TOTAL  AMMONIA  AND  CHLORINE  AT  THE 
GUELPH  WPCP  NECESSARY  TO  MAINTAIN  A  40%  LUZ 

Al lowable  Effluent 
Treatment  Present  Effluent  Concentration  to  Meet 

Level     Concentration    Upstream  Flow  Temperature  Year      PWQO  Target 

m-Vsec    cfs      °C  mg/L 

Tertiary       1.5  mg/L      0.5     17      5 
total  ammonia 

1.8     64      25 


2.5     88      25 


Dechlorination  500  ug/L      0.5      17 


2001M 

1.45 

2031M 

1.20 

2031H 

1.12 

1979 

1.13 

2001M 

0.82 

2031M 

0.56 

2031H 

0.48 

1979 

1.43 

2001M 

1.02 

2031M 

0.66 

2031H 

0.56 

2031M 

338.75 

(ug/1) 

2031H 

235.50 

(ug/1) 

3.2.4  Gait 

A  60%  zone  of  passage  in  the  Grand  River  below  the  Gait  WPCP  is 
attainable  with  regard  to  total  residual  chlorine.  A  chlorine 
concentration  of  500  ug/1  in  the  effluent  does  not  result  in  the 
LUZ  being  extended  to  40%  of  the  total  river  flow,  so  no  plot  of 
the  results  has  been  presented. 

Critical  instream  concentrations  of  ammonia  are  below  the  PWQO 

for  all  the  design  conditions  (Fig.  15).  Results  presented  in 

Table  8  show  that  the  values  of  X   .  for  ammonia  are  15,100 

sceA 

meters  and  greater  for  CAS  treatment  conditions  at  the  Gait  WPCP. 
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TABLE  8.   SHORELINE  DISTANCES  DOWNSTREAM  OF  OUTFALLS  WHERE  THE  PWQO  FOR 
UNIONIZED  AMMONIA  IS  MET 


Streamf 1 ow 

Sewage 

Temperature 

at  Bridgeport 

Distance 

Downstream  of  Ou 

tfall(m) 

Treatment  Plant 

(°C) 

(cfs) 

1979 

200 1M 

2031M 

2031H 

Waterloo 

5 

40 

* 

• 

* 

* 

(CAS) 

150 

10048 

* 

• 

• 

25 

300 

11448 

17548 

* 

* 

400 

9948 

16298 

20448 

* 

500 

8648 

15048 

19348 

* 

600 

7548 

13898 

18248 

20098 

Waterloo 

5 

40 

2082 

4698 

7384 

8548 

(Tertiary) 

150 

70 

222 

477 

1782 

25 

300 

427 

3340 

6148 

7748 

400 

172 

2282 

2848 

4898 

500 

120 

1532 

2690 

3440 

600 

364 

2182 

2490 

Kitchener 

5 

40 

* 

• 

* 

* 

(CAS) 

150 

* 

* 

* 

* 

25 

300 

14250 

* 

* 

* 

400 

13350 

* 

* 

• 

500 

12500 

* 

* 

* 

600 

11650 

• 

* 

* 

Kitchener 

5 

40 

3300 

5000 

6500 

7150 

(Tertiary) 

150 

650 

1350 

2000 

2350 

25 

300 

3050 

6750 

11350 

14100 

400 

1900 

4400 

7700 

9550 

500 

1350 

3100 

5600 

7000 

600 

1050 

2300 

4250 

5350 

Kitchener 

5 

40 

1600 

2250 

2950 

3250 

(Advanced) 

150 

325 

650 

1000 

1100 

25 

300 

1600 

3200 

5500 

6750 

400 

950 

2050 

3700 

4600 

500 

650 

1600 

2600 

3300 

600 

500 

1100 

1950 

2500 

Guelph 

5 

17 

4025 

5100 

7150 

8200 

(Tertiary) 

25 

64 

5050 

4425 

6050 

6900 

88 

4200 

4075 

5800 

6650 

Gait 

5 

40 

39850 

* 

* 

* 

(CAS) 

150 

26750 

43700 

* 

• 

25 

300 

17200 

22850 

26550 

27700 

400 

16600 

22650 

26500 

27700 

500 

15800 

22200 

26300 

27500 

600 

15100 

21750 

26000 

27250 

*  The  predicted  distance,  *sceA  exceeds  the  distance  to  the  next  downstream 
water  pollution  control  plant. 

-  41  - 


3.2.5  Brantford 

As  stated  elsewhere  in  this  report,  the  two  dimensional  stream 
tube  model  was  not  applied  to  the  Brantford  WPCP  as  the  effluent 
mixes  completely  across  the  river  in  a  short  distance  below  the 
outfall.  Therefore,  concentrations  of  un-ionized  ammonia  and 
chlorine  immediately  below  the  outfall  were  predicted  based  on 
the  instantaneous  completely  mixed  assumption  (ICM)  according  to: 


Ca  «  (Ce)(Qe) 


% 


(7) 


C  =    completely  mixed  concentration  of  pollutant  in  river 

a 

just  below  outfal 1 . 

C  =    concentration  of  pollutant  in  effluent, 
e 

Q  =    effluent  flow  rate, 
e 

Q  =    total  river  flow  rate,  given  by  the  sum  of  upstream  and 
effluent  flow  rates. 


Concentrations  of  ammonia  and  chlorine  immediately  below  the 
outfall  of  the  Brantford  WPCP,  under  various  design  conditions, 
are  presented  in  Tables  9  and  10.  For  each  design  case,  the 
downstream  distance,  X  ,  at  which  the  pollutant  concentration 
is  equal  to  the  PWQf),  was  calculated  from  the  relationships: 

Cst  =  Ca  exp  (-K  Xs/uav)  (8) 

Xs=Uavl09e  (  Ca  \  (9) 


Cst 


where 


C  t  =   the  desired  instream  completely  mixed  concentration 

(mg/L) 
C   =   the  instantaneous,  completely  mixed  concentration 

3 

immediately  below  the  outfall  (mg/L) 
X   =   downstream  distance  from  the  outfall  where  the 

concentration  is  C  t(m) 
K   =   decay  rate  coefficient  (sec"  ,  base  e) 
uav  =   average  stream  velocity  (m/s) 
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TABLE  9:   CONCENTRATIONS  OF  TOTAL  RESIDUAL  CHLORINE  BASED  ON  INSTANTANEOUS 

COMPLETELY  MIXED  ASSUMPTION  AND  DISTANCES  WHERE  PWOO  TARGET  IS  MET 
FOR  THE  BRANTFORD  WPCP  WITH  C  =  500  ug/L 


Population 

Season 

Upstream 
Case 

Design* 
Flows 

Instantaneous 
Completely  Mixed 
Concentration  of 
TRC  (ug/L) 

Downst 
Where 

ream  Distance 
PWQO  Target 
is  Met 

n-Vsec 

cfs 

(m) 

1979 

Winter 

13.8 

486.0 

17.0 

11.6 

16.6 

586.0 

14.0 

10.5 

Summer 

20.5 

726.0 

12.0 

9.7 

23.4 

826.0 

10.0 

8.7 

2001M 

Winter 

15.1 

535.0 

22.0 

13.0 

18.0 

635.0 

19.0 

12.2 

Summer 

22.0 

775.0 

15.0 

10.9 

24.8 

875.0 

13.0 

8.7 

2031M 

Winter 

16.9 

596.0 

30.0 

14.6 

19.7 

696.0 

25.0 

13.6 

Summer 

23.7 

836.0 

21.0 

12.7 

26.5 

936.0 

19.0 

12.2 

2031M 

Winter 

17.3 

611.0 

40.0 

16.2 

20.1 

711.0 

31.0 

14.8 

Summer 

24.1 

851.0 

26.0 

13.9 

26.9 

951.0 

23.0 

13.2 

*  Upstream  flows  at  Brantford  are  a  combination  of  flows  from  the  main  Grand, 
Speed  and  Nith  Rivers. 
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TABLE  10:  CONCENTRATIONS  OF  AMMONIA  RASED  ON  ICM  AND  DISTANCES  WHERE  PWQO 
TARGET  IS  MET  FOR  THE  BRANTFORD  WPCP  WITH  C  =  12.5  mg/L. 


Instantaneous 

Downstream  Distai 

Population 

Upstream 

Flow  at 

Completely  Mixed 

Where  PWQO  Targi 

Season 

Bridgeport 

Concentration  of 
TRC  (ug/L) 

is  Met 

m-Vsec 

cfs 

(m) 

1979 

Winter 

13.8 

486.0 

0.43 

991 

16.6 

586.0 

0.36 

685 

Summer 

20.5 

726.0 

0.29 

312 

23.4 

826.0 

0.26 

124 

2001M 

Winter 

15.1 

535.0 

0.55 

1409 

18.0 

635.0 

0.46 

1126 

Summer 

22.0 

775.0 

0.38 

778 

24.8 

875.0 

0.34 

586 

2031M 

Winter 

16.9 

596.0 

0.74 

1927 

19.7 

696.0 

0.64 

1669 

Summer 

23.7 

836.0 

0.53 

1351 

26.5 

936.0 

0.48 

1181 

2031M 

Winter 

17.3 

611.0 

0.89 

2239 

20.1 

711.0 

0.77 

1996 

Summer 

24.1 

851.0 

0.65 

1703 

26.9 

951.0 

0.58 

1507 

*Based  on  a  pH  of  8.2  and  a  temperature  of  5  C  (winter)  and  25  C  (summer), 


7657g 


44 


The  results  based  on  the  instantaneous  completely  nixed 
assumption  indicate  that  the  PWQO  in  the  Grand  River  below  the 
Brantford  WPCP  outfall  would  be  met  within  very  small  distances 
in  the  case  of  TRC  and  distances  ranging  up  to  2239  meters  for 
un- ionized  ammonia. 

It  can  be  noted  from  a  previous  study  on  the  Brantford  WPCP  by 
Wisz  et  a]_   (1978),  the  total  residual  chlorine  in  the  effluent 
averaged  1399  ug/L.  Residual  chlorine  was  observed  at 
approximately  600  m  downstream  of  the  outfall  and  that  some 
lateral  variability  in  the  TRC  concentrations  occurred.  Based 
on  a  comparison  of  this  field  data  with  the  predictions 
presented  in  Table  9,  it  appears  that  the  decay  rate  used  in  the 
instantaneous  complete  mixing  calculations  was  somewhat  high. 
However,  for  the  sake  of  consistency,  the  same  decay  rate  has 
been  used  in  all  five  Grand  River  studies  presented  in  this 
report. 
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3.3  Assumptions  and  Limitations  of  the  Model 

The  model  used  in  this  study  is  applicable  only  to  steady  state 
conditions. 

A  major  assumption  of  the  model  is  that  no  significant 
background  concentration  of  the  pollutant  occurs.  The  model 
only  considers  the  effect  of  pollutant  concentration  in  the 
effluent  released  from  the  wastewater  treatment  plant  under 
consideration. 

Decay  rates  associated  with  total  residual  chlorine  can  be 
highly  variable  from  one  location  to  another.  Previous  studies 
(Gowda,  1980a)  have  found  TRC  decay  rates  in  streams  to  lie  in 
the  range  of  0.56  x  10"4  to  73  x  10"4  per  second  at  20°C 
(to  base  e).  Since  TRC  decay' rates  were  not  measured  at  each 
site  considered  in  this  study,  the  decay  rate  used  in  the  model 
was  an  average  of  these  observed  extreme  rates.  Results  of  a 
sensitivity  analysis  on  the  chlorine  decay  rates  at  the  Waterloo 
WPCP  indicated  that  with  the  low  TRC  decay  rate,  the  PWQO  for 
chlorine  is  violated  in  1979  (300  cfs,  25°C).  With  the  higher 
decay  rate,  no  violations  occurred.  Therefore,  a  more  detailed 
examination  of  chlorine  decay  rates  within  the  basin  should  be 
undertaken. 

As  stated  previously,  the  Kitchener  WPCP  outfall  is  located  at 
the  mid-point  of  the  river.  Since  the  model  is  based  on  a  bank 
outfall,  the  outfall  for  the  Kitchener  plant  has  been  considered 
to  be  moved  to  a  different  location  for  modelling  purposes. 
More  thorough  studies  would  have  to  be  conducted  to  determine  if 
an  improvement  in  effluent  quality  would  alleviate  the  problem 
of  the  effluent  barrier  which  occurs  with  the  current  discharge 
location. 
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Due  to  data  limitations  and  time  and  resource  constraints,  the 
model  was  validated  for  a  conservative  material  for  the  Waterloo 
site  only.  Therefore,  the  results  presented  in  this  report 
should  be  considered  to  be  indicative  only  of  general,  relative 
trends  in  water  quality  changes  resulting  from  various 
management  options. 

More  detailed  descriptions  of  model  assumptions  are  given  in 
Water  Resources  Paper  #14  (Gowda,  1980b). 
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IV  SUMMARY 

Results  of  application  of  the  mixing  zone  model  indicate  that  within 
the  Grand  River  basin,  for  the  design  cases  considered,  a  zone  of 
passage  consisting  of  60%  of  the  total  river  flow  is  attainable  with 
regard  to  chlorine  toxicity.  Only  beyond  the  year  2001  projected 
population  level  at  Guelph,  under  winter  conditions,  are  the  PWQO 
for  chlorine  violated  outside  the  mixing  zone. 

In  order  to  maintain  the  0.02  mg/L  objective  for  un-ionized  ammonia, 
at  the  boundary  of  the  mixing  zone,  tertiary  treatment  is  necessary 
at  the  Waterloo  wastewater  treatment  plant.  With  implementation  of 
tertiary  treatment  at  Kitchener,  Provincial  Objectives  for 
un-ionized  ammonia  can  be  achieved  until  2031  (summer  low  flow);  and 
with  the  subsequent  addition  of  advanced  treatment  (ACA),  no 
violations  of  the  objective  would  occur.  The  current  tertiary 
treatment  level  at  the  Guelph  wastewater  treatment  plant  is  not 
sufficient  to  achieve  the  0.02  mg/L  objective  under  summer  low  flow 
conditions.  At  Gait,  no  additional  treatment  over  conventional 
activated  sludge  is  necessary  to  comply  with  the  PWQO  for  un-ionized 
ammonia. 

Due  to  complete  cross-sectional  mixing  of  effluent  in  the  river  in  a 
short  distance  below  the  outfall,  the  two  dimensional  mixing  zone 
model  was  not  applied  to  the  Brantford  case.  However,  results  based 
on  the  instantaneous  complete  mixing  assumption  indicate  that 
residual  chlorine  concentrations  may  exceed  the  PWQO  for 
approximately  20  m  downstream  and  un-ionized  ammonia  concentrations 
may  exceed  objectives  for  124  to  2239  m  below  the  outfall,  depending 
on  the  design  conditions. 

Given  the  limitations  previously  identified,  the  predicted  results 
should  be  considered  to  be  indications  of  water  quality  changes  due 
to  various  management  options.  Detailed  field  and  modelling  studies 
at  appropriate  implementation  stages  are  recommended,  in  order  to 
update  the  results  presented  in  this  report. 
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A-l         LIST   OF   MIXAPPLN    PROGRAM 


OOO00O10 

00600020 

00000030 

00000040 

00000050 

OOOOOOoO 

00000070 

00000080 

00000090 

00000100 

000001 10 

00000120 

06000130 

00000140 

06060150 

00066160 

00000170 

00000180 

00000190 

00000200 

00000210 

00000220 

60000230 

00000240 

00000250 

00000240 

0060026I 

00000270 

00000280 

00000290 

00000300 

00000310 

00000320 

00000330 

00000340 

00006356 

00000366 

06666370 

00000380 

00000390 

00000400 

00000410 

00000420 

00000430 

00000440 

00000450 

00000460 

00000470 

66006480 

00000490 

00006500 

00000510 

00000520 

•J0000530 

60000540 

00000550 

00000560 

00000570 

00000580 

00000590 

00000600 

00000610 

00000620 

00000630 

00000640 

00000650 

00666660 

00000670 

00000680 

00000690 

60000700 

00000716 

00000720 

00000730 

00000  74  0 

.10000750 

00000760 

00066770 

00000786 

60660790 

00000800 


C      PROGRAM  NAME    MIXAPPLN   »  »  STREAMTUBE  MODEL  FOR  PIPE  OUTFALL 

C     THIS  PROGRAM  IS  SET  UP  FOR  CONSERVATIVE,  NONCONSERVATIVE  UITH 

C     FIRST  ORDER  DECAY ( VIZ ., RESIDUAL  CHLORINE.  PHENOL.  RADIONUCLIDES, 

C     INDICATOR  BACTERIA),  AND  UN-IONIZED  AMMONIA  CONSTITUENTS. 

C     THIS  PROGRAM  INCLUDES  OPTIONS  FOR  DESIGN    QRI  VER ,  QEFFL  .  TEMP  i,  PH. 

C      PROGRAM  DEVELOPED  BY  T.  P.  H.  GOUDA,  UATER  RESOURCES  BRANCH 

C       DATE   JUNE  1980. 

C   DECLARATION  STATEMENTS 

C 

DIMENSION  TITLE<26) , C( 1 0, 50) , CUI < 1 0,50) , PH( 4 ) 

DIMENSION  ARAY1 (400, 10) , ARAY2( 400 , 4 > , ARAY3< 466 , 4 ) 

REAL»8  X(10) ,XX(10) . PBK ,  TBK . PXUC, TXUC , A3 , RKS( 1 0 > . QY(50> , QRS . 
1THETA,BFUR,HPUR,UPUR.QRT0,QT,QRUP(6> ,QEFL(6) , CTDP , PHDG , VOL ( 10) , 
2B<10> ,H( 10)  U(  10) ,BSUM( 10) ,TOT< 10) , TMP(6> , RQ( 1 1  > , TEMPS, PAX  1  . FAX2, 
3RF<  10)  ,BS(  10)  ,HS(  10)  ,US<  10)  ,  BETAdO)  ,  BU(  1  6  )  ,  HU(  1  0  )  ,  UU(  1  0  )  ,  RKT  (  1  6  )  , 
4PHK 10) ,R,QCP.PQX,EY(16) , XCRIT , PCRIT , XEK , XMZ , PMZ , TMZ , RBK . RCR T , REK , 
5RXS  ,  RK  AV  ,  AUCP  ,  RFUC  ,  XUCP  ,  RBT  ,  RBKG  ,  XSCE  ,  XSCEA  ,  PAX1  ,  PAX2  .  PUI1  .  PU2  ,  QCR  . 
6PHUC. DELQ.XL 

CALL  ERRSET(  208,  999,-1  ) 
C 
C   INPUT  DATA 

37  URITE<6,2> 

READd  .3)  TITLE 

URITE(6, 4) 

READd  ,  »)     QRS,BPUR,HPUR,UPUR,  TEMPS,  NTR 

URITE(4,6> 

DO    15    1=1 ,NTR 
15    READd,»)     X<  I  )  ,BS(  I)  ,HS(I)  ,US(I) 

URITE(6,8) 

READd.»)     (BETA(I)  ,  1  =  1  ,NTR) 

WRITE (6, 30) 

READd,»)    MQ,  <QRUP( J) . J»1  ,MQ> 

URITE(6,82> 

READ( 1 ,»)MT, (TMP(L) ,L"1 ,MT) 

URITE(6,86) 

READd,»)    MF,  (QEFL(L)  ,L=1  ,  MF  ) 
35       URITE     (6,52) 

READ( 1 ,»)     AMONIA 

IFCAMONIA.LE.OIGO    TO    70 

URITE (4, 57) 

READ( 1 , «)MPH, (PH( JPH) , JPH=1 ,MPH) 
70   CONTINUE 

URITE(6,56) 

READ (1 . » ) QCP . CEFL , CBKG , CS , THETA , RBK , XUCP 

URITE(6,55) 

READd,»)  (RKS(I)  ,  1  =  1  ,NTR) 
C   CALCULATE  FLOU  &  TEMP'R  SCALE-UP  PARAMETERS 
C     DELTA=0.OOO1  ,    ALOGd/DELTA)  =  (4.0»2.3O26)  FOR  CALCN.  OF  BGX 

NQR=1 1 

NU=6 

DO  10  L=1 .NQR 
10   RQ(L>=(L-1 )/10.0 

URITE(6,24)  TITLE 

IRUN=0 
C      BEGIN  COMPUTATIONS  FOR  THE  INPUT  OPTIONS. 

DO  20  JF=1 ,MF 

DO  20  JQ=1 ,M8 

QT=QRUP ( JQ ) *QEFL  <  JF ) 

DELQ=QT/(NQR-1 ) 

QRTO=QT/QRS 

CA=(CEFL»QEFL< JF)/QT) 

DO  26  JT=1 ,MT 

CTDP=THETA»»(TMP< JT) -TEMPS) 

RBT=RBK»CTDP 

IF (AMONIA. LE. 0)MPH=1 

DO  20  JPH=1 ,MPH 

IRUN=IRUN*1 

IR=IRUN 

IF < AMONIA. LE. 0)PH( JPH) =7.0 

ARAY1 ( IR, 1 )=IRUN 

ARAY1 (IR,2)=QEFL( JF) 

ARAY1 <IR,3)=QRUP( JQ) 

ARAY1 ( IR, 4)=TMP( JT) 

ARAY1 (IR,5)=PH( JPH) 

UMTE(6,92)IRUN 
C      CALCULATE  AMMONIA  IONIZATION  PARAMETER 

IF(AMONIA.LE.O)  GO  TO  34 

PK  A=0  .  0901  8*2729  .  92/  (  TMP  (  JT  )  -»273  .  2  ) 

PF=PKA-PH( JPH) 

PCTU=1  .  /d  .»10.»»PF) 

CSTTL=CS/PCTU 


LIST   OF   MIXAPPLN    PROGRAM     (CONTINUED) 


VALUES  BU.HU.UU  FROM  OUTFALL  TO  TRANSECT(I) 


0OOOO810  34   URITE<6,88>0EFL<  JF>  ,ORUP(  JQ)  ,TMP<  JT>  , PH< JPH > . CEFL . CS 

00090820  IFCAMÛNIA.GE. 1 )URITE< 6 . 54 JCSTTL 

0G0G0B30  URITE<6.90>  t RQ( L ) . L  =  1 , NU ) 

00000840  DO  12  I-t,NTR 

ooeeosso  bsum<i>=o. 

00000860  VOL(t>=0.0 

00000870  12   TOT(I>=0. 

00000880  C  BEGIN  COMPUTATIONS  AT  TRANSECT . I . 

00000890  DO  14  I^I.NTR 

00060900  C  CALCULATE  B,H,U  FOR  DESIGN  FLOU=BT,  FROM  LEOPOLD-MADDOCK  EONS. 

00000910  B< I>=BSC I >»QRTO»«BPUR 

06000920  H(I)=HS( I >»CRTO«»HPUR 

00000930  U< I)=US< I>»QRTO»»UPUR 

00000940  C  CALCULATE  WEIGHTED  MEAN 

00000950  IF(I.GE.2>  GO  TO  60 

00000960  XXtt)»X(tl 

00000970  BU( 1 >=B< 1 ) 

90000980  HU( 1 >=H< I  ) 

00006990  UU(  I  )=U( 1  ) 

00801000  BSUMd  >=BM  >»XX(1  > 

00001010  VOL< t >=XX< I >»B< 1 >»H(1 > 

00001020  T0TC1 )=XX(1 >/U<1 ) 

00001030  GO  TO  62 

00001040  60  11=1-1 

•10001050  XX<I>=X< I)-X(I1 > 

00001060  BSUM(I)=BSUM<It ) ♦  © . 5»XX< I ) • ( B( I 1 )»B(I) ) 

00001070  V0L<I)=V0LU1 >*0.25«XX< I)«<B(I1 )*B(I>)«<H(I1 >+H<I) ) 

jG0O1 080  TOT( I)=T0T(I1 )+XX<I)/U( I) 

00001090  BU< I >=BSUM(I)/X<I> 

00001  100  HU<  I)=VOL(  I)/(X(  I)»BU(D) 

10001 110  UU( I)=QT/(BU( I )«HU( I) ) 

30001120  C  CALCULA1E  PRODUCT  FUNCTION  FOR  DECAY,  RF ( I ) ,  &  DISPERSION  FACTOR 

30001130  62  RKT<  I  >=CTDP«RKSC  I.) 

30001140  PBK=RBT»TOT< I) 

00001150  CALL  PDET(PBK , TBK 1 

00001160  CBKX=CBKG«TBK 

06601  '.  ?0  A3=IRKT< I >»XX(I) )/U(I> 

00001180  R=D£XP<-A3> 

00001190  IFCI.GE.2)  GO  TO  64 

00601 JOO  RFC  1  >=R 

00061210  GO  TO  66 

90061220  64  RF( I )=RF( 11 )»R 

00001230  66  EY( I)=BFTA( I)»B( I)«U(I) 

00001  L?40  PHI(I)=BETA<  I)»X(  I)/BU<  I) 

0006125ft  PHDG=4.0«PHI(I) 

•J0O0 1  260  CMAX=0  .  5«CA»RF  (  I  )  /DSQRT  (3.141  6«PH  I  <  I  )  ) 

00001270  C  LATERAL  CONC.  DISTR'N  AT  TRANSECT , I . 

•0001280  DO  16  K=1 ,NU 

00001290  QY<K >=<K-1 )»DELQ 

00001300  IF  (QY(K).CT.BT)  QYCK)=QT 

06061310  PAX1=(QY(K1-QCP)/QT 

00001320  PAX2=(QY(K)tQCP>/QT 

06001330  CALL  SUMSRSiPAXI . PAX2 , PHDG , SUMT ) 

00061340  C<I .K >=CMAX»SUMT*CBKX 

00001350  16   CONTINUE 

00601360  URITE<6,25>  X< I > . EY ( I ) , ( C( I . K > , K  =  1 , NU ) 

00001370  14   CONTINUE 

00001380  C  UN-IONIZED  AMMONIA  CONCENTRATION  DISTRN. 

00601390  IF(AM0NIA.LE.6)G0  TO  19 

O60O146G  URITEC6.26) 

O0001410  DO  19  1=1 ,NTR 

•36601420  DO  17  K  =  1,NU 

06601430  17   CUKI,K)=C(  I,K)»PCTU 

00001440  URITE16.25)  X(  I  >  ,  EY<  I  >  ,  <  CUI  <  I  ,  K  >  ,  K>=1  .  NU  > 

60061450  19   CONTINUE 

00661460  C  END  OF  COMPUTATIONS  AT  TRANSECT , I . 

60601470  IF(AMONIA.LE.0)CSL=CS 

■06061480  IFtAMONIA.GE.  1  )CSL=CSTTL 

00061490  C  COMPUTE  BACKGROUND  AVG .  CONC.  AT  D/S  UPCP . 

00001500  RKAV=-DLOG(RF(NTR) )/TOT<NTR) 

106011.10  AUCP=RKAV»XUCP/UU(NTR1 

00001520  CALL  PDET'AUCP.RFUC) 

O60O153O  CBA=CA»RFUC 

00601540  76   PXUC=RBT«XUCP/UU(NTR) 

30661556  CALL  PDET (PXUC . TXUC ) 

06001560  CBB=CBKG»TXUC 

000015  70  CAl'P=CBA*CBB 

00001580  C  COMPU.~  BANK  CONC.  AT  D/S  UPCP. 

06061590  PU1=0. 

00001600  PU2=0. 

00061610  PHUC=4.0«BETA(NTR)»XUCP/BU(NTR) 


LIST    OF   MIXAPPLN    PROGRAM     (CONTINUED] 


00001620 

00001630 

00001640 

00001650 

00001660 

00001670 

00001680 

00001690 

00001700 

00001710 

00001720 

00001730 

00001740 

00001750 

00001760 

00001770 

00001780 

00001790 

00001800 

00001810 

00001820 

00001830 

00001840 

00001850 

00001860 

00001870 

00001880 

00001890 

00001900 

00001910 

00001920 

00001930 

00001940 

00001950 

00001960 

00001970 

00001980 

00001990 

00002000 

00002010 

00002020 

00002030 

00002040 

00002050 

00002060 

00002070 

00002080 

00002090 

00002100 

000021 10 

00002120 

00002130 

00002140 

00002150 

00002160 

00002170 

00002180 

00002190 

00002200 

00002210 

00002220 

00002230 

00002240 

00002250 

00002260 

00002270 

00002280 

00002290 

00002300 

00002310 

00002320 

00002330 

00002340 

00002350 

00002360 

00002370 

00002380 

00002390 

ooco;  ^o 

000024.  i 
00002420 


CALL  SUMSRSIPW1  , PW2 , PHUC , SUMUC > 

CWCP=CBA»SUMUC/DSQRT(3. I  41 6»PHUC ) *CBB 
:   COMPUTATIONS  FOR  MIXING  ZONE  PARAMETERS 

XMZ=B< NTR) /BETA (NTR) 

PMZ=RKAV»XMZ/UU(NTR> 

CALL  PDETCPMZ.TMZ) 

CMZ-CA»TMZ 
2      CALCULATE  XSCE 

IF(CUCP.GT.CSL)  GOTO  220 

CALL  PARSPR ( C . CSL , X , NTR , KXS , XEST , CXS ) 

RXS=-DLOG(RF(KXS> ) /TOT (KXS) 

CALL  SPREAD ( CSL , CA , CBKG , RBT , CXS , XEST , XSCE , KXS , BU , RXS . UU , BETA , ITRN ) 

GO  TO  222 
220  XSCE=-999.0 

ITRN=0 
222  URITE(6,50)XSCE, ITRN , XMZ , CMZ . XUCP , CUCP , CAUP 

ARAY1 < IR,6)=CAUP 

ARAY1 ( IR.7)=XSCE 
C   CRITICAL  POINT  RESULTS  FOR  QRATI0=0.2  ,  0.3  4  0.4 

URITE(6,94) 

DO  20  K-3,5 

KK=K-2 

CCRIT=-1 .OE+10 
C   SEARCH  FOR  TRANSECT  NEAR  WHICH  CRIT.  CONC.  OCCURS  TO  FIND  MOVING 
C      AVG.  VALUES  FOR  CRITICAL  POINT  COMPUTATIONS 

DO  18  1=1 ,NTR 

IF(CCRIT.GE.C(I,K))  GO  TO  18 

ICR=I 

CCRIT=C(I,K> 
18   CONTINUE 

XCRIT-=X(  ICR) 

PCRIT-PHI  '.  ICR) 

RCRT=-DLOG(RF(ICR))/TOT(ICR> 
C    CALCULATE  XL.CL  (.  CEA  BY  CRIT.  POINT  METHOD  USING  MOVING  AVG.  VALUES 
40   QCR=4.0»RCRr»XCRIT/(PCRIT«UU< ICR) > « < QY ( K  > /QT >«»2 

XL=0.25»UU(ICR>»('-1  .O+DSQRTO-fQCR)  )/RCRT 

XLDIF=DABS'XL-XCRIT) 

XPCT=100«XLDIF/XL 

IF<XPCT.LE.5. )  GO  TO  42 

XCR  I  T  =  XL 

PCRIT=BETA(ICR)»XL/BU(ICR> 

GO  TO  40 
42   PQX=><RCRT«XL/UU(ICR>>  +  (0.25/PCRIT»<QY(K>/QT>»«2> 

CL=(CA/DSQRT<3. 1416«PCRIT> )»DEXP( -PQX) ♦CBKG»DEXP< -RBT»XL/UU( ICR) > 

CEA=CEFL»CSL/CL 
C      COMPUTE  XSCEA 

IF(CUCP.GT.CL)GO  TO  114 

CALL  PARSPR  <  C , CL . X . NTR , IEK . XEK . CXK ) 

REK=-DLOG(RF(IEK)>/TOT(IEK> 

CALL  SPREAD  <  CL . CA , CBKG , RBT , CXK , XEK , XSCEA , IEK , BU , REK , UU , BETA , I TRA ) 

GO  TO  116 
114  XSCEA=-999.0 

ITRA=0 
116  URITE(6.96)RQ<K) . XL , CL , CEA, XSCEA, ITRA 
C      STORE  OUTPUT  IN  ARRAYS 

GO  TO  (45,46,47) ,KK 

45  ARAY1 < IR,8)=RQ(3> 
ARAY1 <IR,9)=CEA 
ARAY1 < IR, 10)»XSCEA 
GO  TO  20 

46  ARAY2< IR, 1 )=RQ(4> 
ARAir2<  IR,2)=CEA 
ARAY2( IR,3)=XSCEA 
GO  TO  20 

47  ARAY3UR,  1  >=RQ(5> 
ARAY3(IR,2)=CEA 
ARAY3< IR,3)=XSCEA 

20   CONTINUE 
C      END  OF  COMPUTATIONS  FOR  THE  INPUT  OPTIONS 
C     PRINT  SUMMARIZED  RESULTS 

URITE(6,78)TITLE 

URITE(6, 100) 

DO  48  N=t , IR 

WRITE (4, 102) (ARAY1 (N, J) , J»1 . 1 0> 

WRITE (6. 104) (ARAY21N, J) , J=1 ,3) 

48  WRITE (6, 104) (ARAY3CN, J) , J=1 .3) 
C   CHANGE  PARAMETERS 

WRITE(6,2B) 
READd  ,»)  NCHNG 
GO  TO (35. 37, 999) , NCHNG 
C   FORMAT  STATEMENTS 


LIST    OF   MIXAPPLN    PROGRAM     (CONTINUED) 


000024  30 
00002-410 
00002450 
00002460 
00002470 
000024BO 
00002490 
00002500 
00002510 
00002520 
00002530 
00002540 
00002550 
00002560 
00002570 
Ô00O25BO 
00002590 
00002600 
00002610 
00002620 
00002*30 
00002640 
00002650 
00002660 
00002670 
00002680 
00002690 
00002700 
00002(10 
00002720 
00002730 
90002740 
00002750 
00002760 
00002770 
00002780 
00002790 
00002800 
00002810 
00002820 
00002830 
00002840 
00002850 
00002860 
00002870 
00002880 
00002990 
00002900 
00002910 
00002920 
00002930 
00002940 
00002950 
00002960 
00002970 
00002980 
00002990 
00003000 
00003010 
00003020 
00003030 
00003040 
00003050 
00003060 
000030  70 
00003080 
000O309O 
00003100 
000031 10 
00003120 
00003130 
00001140 
00003150 
00003160 
00003171 
000031 80 
00003190 
00003200 
00003210 
00003220 
00003230 


3 

4 
6 
8 

24 
25 
26 
28 

50 


52 
54 
55 
56 
57 
78 
80 
82 
86 
88 

90 
92 
94 

96 
100 


104 
99r> 


FORMAT! /'  ENTER  TITLE  OF  STUDY') 

FORMAT*  20A4) 

FORMAT  ('  ENTER  ORS  ,  BPUR  ,  HPUR  ,  UPUR  ,  TEMPS  ,  NTR  '  ) 

FORMAT !'  ENTER  NTR  VALUES  OF  X , BS , US . US • ) 

FURMAT<'  ENTER  NTR  VALUES  OF  BETA') 

F0RMATI2A4» 

FORMAT! 1H1//6X, 'PREDICTIONS  OF  RUNS  FOR  MANAGEMENT  OPTIONS ' /20A4 ) 

F0RMAT(2X,F8. 1 ,1X,F7.4,11F9.3) 

F0RMAT(2X, ' TOXIC  AMMONIA   ) 

FORMAT!/'   TO  CHANGE  PARAMETERS  ENTER  THE  NUMBER  STATED  BELOW/ 
*'    AMMONIA   1,    TITLE   2,    END   3.') 

F0RMAT</5X, ' XS  i  WITH  CE  )  =  '  ,  F8 . 1  ,  3X  ,  '  NO  .  OF  ITERATIONS»' , 12, /5X, 
•'MIXING  ZONE  LENGTH=" ,F9. 1 ,5X, 'CONC=' .F8.2/5X, 'DIST.  TO  D/S  UPCP= ' 
•,F8. 1 ,5X, 'SHORE  CONC.  AT  D/S  UPCP= ' , F6.2/8X , ' AVG .  CONC.AT  D/S  UPCP 
•=  ' ,F6.2) 

FORMAT!'    UN-IONIZED  AMMONIA:  ENTER  1  FOR  YES;  0  FOR  NO') 

F0RMAT(5X,'   CRITERION  FOR  TOTAL  AMMONIA,  CS=',F6.3> 

F0RMAT(2X,'  ENTER  NTR  VALUES  OF  RKS') 

FORMAT!'    ENTER  QCP , CEFL , CBKG , CS , THETA , RKB , XUCP ' ) 

FORMAT!'    ENTER  MPH  AND  PH  VALUES') 

FORMAT! 1H1//6X, 'SUMMARY  OF  RUNS  FOR  MANAGEMENT  OPTIONS ' /20A4 > 

FORMAT!'   ENTER  MQ  J.  QRUP  VALUES') 

FORMAT!  '   ENTER  MT  J,  TMP  VALUES') 

FORMAT!  '   ENTER  MF  J,  QEFL  VALUES') 

F0RMAT(2X,  'I3EFL='  ,F7.3,2X,  '  QRUP=  '  ,F7.3,2X,  '  TEMPR- ' 
»F4.1 ,2X, 'CEFL=- , F7 . 2 , 2X • CS= ' ,F5.2) 

FORMAT  I /6X,  'X' ,6X,  ' EY '  ,7X,6(F3. 1  , 6X ) , / ) 

FORMAT!/'    •  •  RUN  NO.   ',14) 

F0RMATI/9X,'     CRITICAL  POINT  METHOD  RESULTS 
•  '  XL  '  ,  8X  ,  '  CL  '  ,  9X  ,  '  CEA  '  ,  6X  ,  '  XSCEA  '  ,  6X  ,  '  I  TRN  '  ) 

FORMAT (VX.F4. 2, 4X,F6.  1  ,  3X  ,  F7  .  3  ,  4X  ,  F8  .  2  ,  3X  ,  FB  .  2  ,  4X  ,  14  ) 

FORMAT !//1 / . 'RUN*' .3X, 'QEFL' ,5X. ' QPUP ' , 3X , ' TEMP • , 4X , ' PH ' , 4X . 
1 'CAUP' ,5X, 'XSCE' ,4X, 'QY/QT' ,4X, 'CEA' , 5X , 'XSCEA'/) 

FORMAT!  1X,F4.0,2X,,F6.3,2X,F6.2,2C2X.F5.1  )  ,2X,F6.3, 
•2X.F8. 1 ,2X,F4.2,2X,F7.2,2X,F8.1 ) 

F0RMAT!55<,F4.2,2X,F7.2,2X,F8.1 ) 

STOP 

END 


, F4. 1 ,2X, 'PH= 


'/9X, 'QY/fJT' ,6X, 


SUBROUTINE  FOR  SUMMATION  OF  EXPONENTIAL  SERIES  TERMS. 

SUBROUTINE  SUMSRSIQAX1  ,  QAX2  ,  PFDR,  SUM  ) 

REAL»B  QAX1 , QAX2 , PFDR , BGX , SBG , P1 , P2 , P3 , P4 , T1 , T2 , T3 , T4 

BGX=2.3026«PFDR 

SBG=DSQRT(BGX) 
DETERMINE  NO.  OF  IMAGES  REQUIRED 

AN1=(0.5«QAX1 -SBG) -0.5 

AN2=!0.5»QAX1+SBG>+0.5 

AN3=-AN2 

AN4=-AN1 

NM1=IFIX!AN1 ) 

NM2=IFIX!AN2) 

NM3=IFIX!AN3> 

NM4=IFIX(AN4) 

NN1 =1 ♦NM2+IABS1NM1 ) 

NN2=1*IABS!NM3)+IABS(NM4) 

IF!NN1 .GE.NN2)NN=NN1 +1 

IF1NN1 .LT.NN2)NN=NN2+1 
COMPUTE  SUM  OF  EXPONENTIAL  SERIES  TERMS 

SUM=0. 

DO  32  J=1 ,NN 

N-J-1 

P1=!0AX1-2.»N)»»2/PFDR 

P2=(QAX2*2.»N)»«2/PFDR 

CALL  FDET(P1 , T1 ) 

CALL  PDETIP2.T2) 

IFiN.LE.O)  GO  TO  30 

P3=IQAX1 +2. »N)»«2/PFDR 

P4=!DAX2-2.»N)««2/PFDR 

CALL  PDET!P3,T3) 

CALL  PDET(P4,T4) 

GO  TO  32 
30   T3=0. 

T4  =  0. 
32   SUM=SUM*T1*T2*T3*T4 

RETURN 

END 


COMPUTE  EXPONENTIAL  TERMS 
THAT  EXP!-P)=4.3E-18. 
SUBROUTINE  PDEI !P, T) 
REAL'S  P, T 


SET  (P. LE. 40.0)  TO  AVOID  ERROR  208.  SO 


LIST    OF    MIXAPPLN    PROGRAM     (CONTINUED) 


000031' '10 

00003256 

00003260 

00003270 

00003260 

00003290 

00003300 

00003310  C 

00003320  C 

00003330 

00003340 

00003350 

00003360 

00003370 

00003380 

00003390 

00003400 

00003410 

00003420 

00003430 

00003440 

00003450 

00003460 

00003470 

00003480 

00003490  C 

00003500 

00003510  C 

00003520  C 

00003530 

00003540 

00003550 

00003560 

00003570 

00003580 

00003590 

00003600 

00003610 

00003620 

00003630 

00003640 

00OU3A5O 

00003660 

00003670 

00003680 

00003690 

00003700 

00003710 

00003720 

00003730 

00003740 

00003750 

00003760 

00003770 

000037AO 

00003790 

00003800 


IFCP.GE.40. )G0  TO  10 

T-DEXP(-f') 

GO  TO  12 
10   T  =  0.0 
12   CONTINUE 

RETURN 

END 

COMPUTE  PARAMETERS  IN  THE  SUBROUTINE  'SPREAD'. 
SUBROUTINE  PARSPR(C,CL,X,NTR, IEK , XEK ,CXK) 
DIMENSION  C( 1 , 1 ) 
REAL»8   X(1 ) ,XEK 
IF<C<NTR,  1  )  .GE.CL)  GO  TO  75 
DO  72  1=2, NTR 
11=1-1 

IF(C(  1  .  1  )  .LE.CDGO  TO  73 
IF<C(  II ,  1  > .GT.CL.AND.CC I , 1 > .LE.CL)  GO  TO  73 

72  CONTINUE 

73  IEK=II 
GO  TO  79 

75   IEK=NTR 
79   XEK=X(IEK) 

CXK=C( IEK . 1 ) 

RETURN 

END 

SUBROUTINE  SPREAD < CST , CAV , CBG , RB , CXY , XEY , XST , M , B I ,RS,US,BTA, IT) 
COMPUTATION  OF  MAX.  LONGL .  SPREAD, XS,  ALONG  OUTFALL  BANK  UHERE 
(C(XS,0)-CS)=5  PERCENT(ABSOLUTE) . 

REAL»8  Bid  ),BTA<  I  ),RS,UJ(1  )  ,  PS  ,  TS  ,  PHB  ,  RFS  .  RB  ,  XS  ,  PX  1  ,PX2.TG,PG 

CXX=CXY 

XS=XEY 

IT  =  0 
18   IT=IT*1 

DIFF=(CXX-CST) 

RCAB=ABS<DIFF/CST) 

PRCNT=100.»RCAB 

IF<PRCNT.LE.5. )G0  TO  26 

IF< IT.GT.30)  GO  TO  27 

FRX=RCAB/(  1  .O-fPS) 

IF(DIFF.LE.O.O)  XS=XS»( 1-FRX) 

IF(DIFF.GT.O.O)  XS=XS*( 1+FRX) 

PHB=4.0«BTA<M)»XS/BKM> 

PS=RS»XS/US(M) 

CALL  PDET(PS.TS) 

CMM=CAV«TS/DSQRT<3. 1416»PHB> 

PX1=0. 

PX2=0. 

PG=RB«XS/US(M) 

CALL  PDET(PG,TG) 

CALL  SUMSRS(PX1 , PX2 . PHB , SUMP ) 

CXX=CMM»SUMP+CBG»TG 

GO  TO  18 
27   XS=-888.0 
26   XST=XS 

RETURN 

END 


A-2        TYPICAL   INPUT/OUTPUT  OF   MIXAPPLN   PROGRAM 
INPUT: 


ENTER    TITLE    OF    STUDY 

ENTE  R    0  P:  S  ,  B  P  W  ft  ,  II P  U  R  .  U  P  14  ft  .  T  E  M  P  S  .  N  T  ft 
ENTEft    NTft    VALUES    OF    X , BS . US . US 
ENTER    NTft    VALUES    OF    BETA 
ENTER    MO    &    Oft  UP    VALUES 
ENTEft    MT    &    THP    VALUES 
ENTEft    MF    •?,    C'EFL    VALUES 
.    UN-IONIZED    AMMONIA:     ENTEft    i     FOP    YES;     0    FOP    NO 
ENTER    MPI-!    AND    PU    VALUES 
ENTER    fJCP.CEFL,CBKG.CS.  THETA  .  RKB.XWCP 
ENTER    NTP    VALUES    G r    RKS 

TO    CHANGE    PARAMETEPS    ENTER' THE    NUMBER    STATED    BELOW 
AMMONIA       1;  TITLE:     2;  END:     3. 


AMMONIA    IN    MIXING    ZONE-GRAND    RIVER    BELOW    WATERLOO    WPCP    * 
12.54     .05    .5     .45    23    5 
122    54.86    .62    .37 
427    86.87     .62     .63 
1082    44.2    .51     .55 
2240    59.44     .4?     .43 
3398    50.29     .69     .37 
.0025     .0020    .0013     .0008    .0010 

1  1.1 32 

2  20  25 
1  0.307 
1 

1     8.3 

0    12.5    0    0.02    1.106    0    20487 

.0000231     .0000231      ,0000231     .0000231     .0000231 

3 


OUTPUT    OF   MIXAPPLN    PROGRAM 


PREDICTIONS  OF  RUNS  FOR  MANAGEMENT  OPTIONS 
AMMONIA  IN  MIXING  ZONE-GRAND  RIVER  BELOU  WATERLOO  UPCP 


«  •  RUN  NO.       1 
QEFL=   ©.397   QRUP=   1.132   TEMPR=20.0   PH=  9.3 
CRITERION  FOR  TOTAL  AMMONIA,  CS  =  0.272 


CEFL=   12.50  CS=    0.02 


EY 


0.0 


0.1 


0.2 


0.3 


0.5 


122.0 

0.0172 

18. 

832 

12. 

579 

3. 

749 

0, 

498 

0. 

.030 

0 

,001 

427.0 

0.0371 

12 

,  104 

10 

169 

6 

,031 

i 

.525 

0 

.746 

0 

,156 

1082.0 

0.0107 

8. 

.906 

8. 

018 

5. 

851 

3, 

.461 

i 

.660 

0 

.645 

2240.0 

6.0069 

6 

,590 

6. 

,  124 

4. 

,914 

3. 

.405 

*> 

.038 

1 

.053 

3398.0 

0.0063 

4 

.108 

3. 

956 

3 

531 

1 

.923 

2 

.243 

1 

,596 

IX  IC  AMMONIA: 

122.0 

0.0172 

1 

.383 

0 

.924 

0 

.275 

0 

.037 

0 

.002 

0 

.000 

427.0 

0.0371 

0 

,889 

0 

,747 

0 

.443 

0 

.185 

0 

,055 

0 

.011 

1082.0 

0.0107 

0 

.654 

0 

,589 

0 

,430 

0 

.254 

0 

.122 

0 

.047 

2240.0 

0.0069 

0 

.484 

0 

,450 

0 

.361 

© 

.250 

0 

.150 

0 

.077 

3398.0 

0.0063 

0 

.302 

0 

.291 

0 

.259 

0 

.215 

0 

.165 

0 

.117 

XS  (UITH  CE>  =  20259.4    NO.  OF  ITERATIONS»  8 
MIXING  ZONE  LENGTH»   45130.5      CONC=     0.01 
DIST.  TO  D/S  UPCP-  20487.0      SHORE  CONC .  AT  D/S  UPCP= 
AVG.  CONC. AT  D/S  UfCP=    0.25 

CRITICAL  POINT  METHOD  RESULTS 

QY/QT       XL  CL          CEA  XSCEA 

0.20      516.0  5.946         ©.57  2312.38 

0.30     1455.0  3.40©         1.00  3956.84 

0.40     2575.7  2.243         1.52  5735.51 


•  •  RUN  NO.       2 
QEFL  =   0.307   QRUP=   1.132   TEMPR=25.0   PH= 
CRITERION  FOR  TOTAL  AMMONIA.  CS=    ©.196 


8.3   CEFL= 


12.50  CS=    0.02 


EY 


0.0 


0.1 


0.3 


0.4 


9.5 


122.0 

0.9172 

18 

,649 

12 

.457 

3. 

,712 

9 

,494 

0 

.929 

0 

001 

427.0 

0.0371 

1  1 

.815 

9 

.927 

5 

887 

2 

.465 

0 

.728 

0 

,152 

1082.0 

0.0197 

8. 

,392 

7 

.556 

5. 

,514 

3 

.262 

1 

.S64 

0 

608 

224©. © 

0.0069 

5 

.734 

5 

.328 

4 

.276 

-t 

.963 

1 

.773 

0 

.916 

3398.9 

0.0063 

3 

,258 

3 

.137 

"» 

809 

2. 

.318 

1 

.778 

1  , 

,266 

TOXIC  AMMONIA 

122.0 

0.0172 

1 

,902 

1 

.270 

9, 

,379 

9 

.05© 

0 

.993 

0 

,000 

427.0 

9.9371 

1 

,205 

1 

.012 

9. 

600 

0 

.251 

0. 

974 

0. 

.015 

1082.9 

9.9107 

© 

,856 

0 

.770 

e. 

,562 

0 

.333 

9 

,159 

0 

.062 

2240.0 

0.0069 

© 

,585 

© 

.543 

0 

,436 

0. 

,302 

9 

.181 

0 

,093 

3398.9 

0.0©63 

0 

.332 

0 

.320 

9 

.286 

0 

,236 

9 

.181 

© 

.129 

XS  (UITH  CE>  =  14413.0    NO.  OF  ITERATIONS=  9 
MIXING  ZONE  LENGTH»   45130.5      CONC  =     0.00 

DIST.  TO  D/S  UPCP=  20487.0      SHORE  CONC.  AT  D/S  UPCP=   0.05 
AVG.  CONC. AT  D/S  UPCP=    0.05 


CRITICAL  POINT  METHOD  RESULTS 
QY/QT       XL         CL  CEA 

0.2©  481.3  5.665  0.43 
0.39  1284.3  3.092  9.82 
9.49     2221.9      1.872         1.31 


XSCEA 

ITRN 

2240.00 

1 

3474.97 

3 

4907.67 

2 

SUMMARY  OF  RUNS  FOR  MANAGEMENT  OPTIONS 
AMMONIA  IN  MIXING  ZONE-GRAND  RIVER  BELOU  UATERLOO  UPCP  • 


RUN*    QEFL       QPUP    TEMP     PH     CAUP      XSCE     QY/QT 


XSCEA 


1.  0.307     1.13    20.0     8.3    9.248    20259.4   0.20 

0.30 
0.40 

2.  0.30?  1.13         25.0  8.3         0.052         14413.0      0.20 

0.30 
0.40 


0 

.57 

2312. 

.4 

1 

.0© 

3956 

.8 

1 

.  52 

5735 

.  *" 

0 

,43 

224© 

.0 

0 

.82 

3475 

.0 

1 

.31 

4907 

.7 

TO  CHANI.F  PARAMETERS  ENTER  THE  NUMBER  STATED  BELOU 
AMMONIA   1 .     TITLE   2,    END   3. 
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A-4  Decay  Rates  for  Chlorine  and  Ammonia  below  the  major  wastewater  treatment  plants 
in  the  mixing  zone  study 


Chlorine 

Ammonia 
CAS  Treatment 

Ammonia 
Tertiary  Treatment 

Ammonia 
ACA  Treatment 

Waterloo 

3.7  x  lO-3 

2.31  x  lO"5 

2.31  x  10"6 

- 

Kitchener 

3.7  x  lO"3 

3.01  x  lO"5 

3.01  x  10-6 

3.01  x  10-6 

Gait 

3.7  x  10-3 

3.01  x  10-5 

3.01  x  10-6 

- 

Guelph 

3.7  x  10-3 

- 

1.16  x  10-5 

- 

Rrantford 

3.7  x  10-3 

1.16  x  lO"5 

- 

- 

*  per  second  at  20°  C  (to  base  e) 
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